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Abstract 

Experimental  Study  of  Sound  Waves  in  Sandy  Sediment 
Michael  W.  Yargus 

Chair  of  the  Supervisory  Committee: 

Professor  Darrell  R.  Jackson 
Electrical  Engineering 

This  dissertation  describes  experiments  intended  to  help  understand  the  physics 
of  sound  (compressional  waves)  propagating  through  sandy  sediments  (unconsolidated 
porous  media).  The  theory  (using  a  lumped  parameter  model)  and  measurements  (using 
a  reflection  ratio  technique)  includes  derivations  and  measurements  of  acoustic  imped¬ 
ances,  effective  densities,  wave  speeds  (phase  velocities),  effective  pressures,  mode 
shapes,  pressure  reflection  coefficients,  and  material  moduli.  The  results  show  the 
acoustic  impedance  divided  by  the  phase  velocity,  rendering  an  “effective  density,”  is 
less  than  the  total  density  of  the  sediment  (effective  density  =  89%  +  3%  of  total).  The 
results  also  show  the  fluid  in  the  sediment  oscillates  back-and-forth  2.2  +  0.4  times  far¬ 
ther  than  the  sand  in  the  sediment  (mode  shape)  during  the  passing  of  a  sound  wave. 
These  facts  suggest  the  existence  of  Biot  waves  (two  compressional  waves)  in  water- 
saturated  sand. 
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Chapter  1: 

INTRODUCTION 

This  dissertation  implements  theory  and  procedures  designed  to  measure  physical 
parameters  of  acoustic  waves  (also  called  sound  waves,  compressional  waves,  bulk  com- 
pressional  waves,  and  dilatational  waves)  in  sandy  sediments.  This  experimental  study 
was  done  to  help  understand  the  mechanics  of  sound  propagating  through  water-saturated 
sand  and  in  particular  to  look  for  effects  only  seen  in  porous  media.  These  effects  in¬ 
clude  different  particle  oscillation  amplitudes  of  the  sand  and  fluid  (pore  fluid  in  the 
sediment),  a  sound  wave  not  “using”  the  total  mass  of  the  sediment,  and  the  existence  of 
two  compressional  waves  with  different  speeds.  These  are  characteristics  of  the  two  (fast 
and  slow)  Biot  waves. 

Sound  wave  measurements  reported  in  this  dissertation  are  obtained  through  a  re¬ 
flection  ratio  technique,  a  new  technique  that  uses  a  reference  surface  and  a  pressure  re¬ 
lease  surface  to  measure  coherently  averaged  pressure  reflection  coefficients  from  a  sur¬ 
face  of  interest.  In  the  first  of  two  parts,  a  signal  of  interest  (reflected  from  the  surface  of 
interest)  is  divided  by  a  reflected  reference  signal.  The  reference  signal  is  a  measure  the 
sound  wave  being  transmitted,  and  by  dividing  the  reference  signal  into  the  signal  of  in¬ 
terest  a  transfer  function  is  obtained.  This  has  the  advantage  of  setting  a  nearly  perfect 
time  reference  so  the  reflections  can  be  averaged  coherently.  In  addition,  the  two  signals 
from  the  transmitter  do  not  have  to  be  identical  (the  ringing  of  the  transmitting  hydro¬ 
phone  was  actively  canceled  and  the  signal  periodically  needed  to  be  “reshaped”).  In  the 
second  part,  the  transfer  function  from  the  surface  of  interest  is  divided  by  a  transfer 
function  from  a  pressure  release  surface.  The  pressure  release  surface  is  set  up  so  that  the 
ratio  of  the  two  transfer  functions  leaves  only  the  pressure  reflection  coefficient.  This 
technique  has  the  advantage  of  measuring  a  pressure  reflection  coefficient  without  need¬ 
ing  to  know  the  acoustic  properties  (dissipation,  damping,  spreading  loss,  etc.)  of  the  me¬ 
dia  the  sound  travels  through  and  without  needing  to  calibrate  the  transducers  (receiver 
and  transmitter).  With  the  reflection  ratio  technique,  coherently  averaged  pressure  re¬ 
flection  coefficients  are  measured  without  calibrating  the  transducers  and  without  use  of 
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geometric  correction  factors.  By  knowing  pressure  reflection  coefficients,  acoustic  im¬ 
pedances  of  the  sand  can  be  obtained. 

The  mathematical  theory  is  based  on  a  lumped  parameter  model,  a  new  model 
that  visibly  shows  the  physics  of  a  Biot  medium  and  facilitates  use  of  free-body  diagrams 
to  find  the  distribution  of  pressures  and  the  relative  velocities  between  particles  and 
acoustic  waves,  e.g.,  sand  particle  velocity  to  fast  wave  “effective”  velocity.  By  knowing 
the  fast  and  slow  wave  “effective”  velocities,  rigorous  derivations  of  the  fast  and  slow 
wave  acoustic  impedances  are  obtained  for  the  first  time. 

1.1  OUTLINE  OF  DISSERTATION 

This  chapter  includes:  (1)  A  background  of  why  knowledge  of  the  physics  of 
sound  waves  propagating  through  sand  is  of  practical  importance.  (2)  A  short  history  of 
sound  propagating  through  saturated  porous  materials.  (3)  Previous  observations  of  Biot 
waves.  (4)  A  look  at  two  contemporary  controversies.  One  is  whether  slow  waves  exist 
in  unconsolidated  porous  materials,  i.e.,  sediments,  and  another  is  whether  the  claimed 
compression  waves  in  sediments  at  1200  m/s  are  slow  Biot  waves  or  artifacts  from  scat¬ 
tering.  (5)  A  list  of  key  parameters  and  assumptions  that  play  a  role  in  the  dynamic  re¬ 
sponse  of  sediments.  (6)  A  list  of  the  new  contributions  published  in  this  dissertation. 

The  second  chapter  of  this  dissertation  goes  through  the  mathematical  theory  and 
discusses  a  lumped  parameter  model.  The  third  chapter  discusses  the  acoustic  measure¬ 
ments  and  the  reflection  ratio  technique.  The  fourth  chapter  discusses  the  analysis  and 
results;  the  last  chapter  gives  a  discussion  with  conclusions. 

1.2  BACKGROUND 

Understanding  the  physics  of  sound  propagation  into  a  sandy  bottom  of  a  shallow 
water  environment  is  driven  by  a  need  to  improve  techniques  for  finding  and  classifying 
objects  that  are  buried.  Sonar,  traditionally  used  to  locate  objects  in  water  or  map  the 
water-sediment  interface,  can  be  used  to  locate  objects  such  as  mines,  cables,  pipelines 
and  containers  buried  in  sediment.  The  principle  reason  acoustic  wave  propagation  and 
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scattering  have  been  widely  used  to  study  proeesses  and  detect  objects  in  the  oeean  is  the 
limited  range  of  optical  and  other  electromagnetic  wave  propagation  underwater.  For  a 
wide  range  of  seientifie  and  engineering  problems  dealing  with  the  oeean  environment, 
acoustie  wave  propagation  and  seattering  is  the  only  remote  sensing  teehnology  avail¬ 
able.  A  buried  target  ean  be  illuminated  with  maximum  intensity  by  a  sonar  signal  nor¬ 
mally  ineident  on  the  sediment  interfaee.  However,  the  time  taken  to  deteet  buried  ob- 
jeets  ean  become  prohibitively  large  if  only  normal  incidence  is  used.  Using  shallow  in- 
eidence  grazing  angles  can  speed  up  deteetion  operations  considerably.  For  sand  sedi¬ 
ments,  the  compressional  sound  speed  in  the  sediment  is  faster  than  in  the  water  and  the 
wave  transmitted  into  the  sediment  is  refraeted  at  a  smaller  grazing  angle  than  in  the  wa¬ 
ter.  At  the  eritieal  angle,  typieally  at  grazing  angles  of  about  20  to  30  degrees,  the  re¬ 
fraeted  eompressional  wave  is  parallel  to  the  interfaee  and  no  longer  propagates  into  the 
sediment.  There  is  significant  interest  in  mechanisms  by  which  a  compressional  wave  in 
the  water  ean  be  eoupled  into  the  sediment  at  angles  below  the  eritieal  angle  with  energy 
sufficient  enough  for  finding  buried  objects  [24,  27,  34].  It  has  been  proposed  that  the 
slow  Biot  wave  may  provide  such  a  mechanism  [34]. 

This  dissertation  explores  normally  ineident  waves  in  an  attempt  to  understand  the 
nature  of  sound  propagating  through  sandy  sediments.  The  results  of  this  work  are  rele¬ 
vant  to  the  general  problem  including  incidence  at  angles  smaller  than  the  eritieal  angle. 

1.3  HISTORY 

A  complete  theory  for  aeoustie  wave  propagation  in  porous  media  was  developed 
by  Biot  [3,  4,  5]  and  presented  in  a  series  of  papers.  Biot  has  developed  a  comprehensive 
theory  for  the  static  and  dynamic  response  of  linear,  porous  materials  eontaining  com¬ 
pressible  fluid.  This  theory  was  derived  by  adding  inertia  terms  into  the  equation  of  eon- 
solidation  developed  in  earlier  theories.  Biot’s  theory  predicts  that  in  the  absenee  of 
boundaries,  three  kinds  of  body  waves,  two  dilatational  and  one  shear,  may  exist  in  a 
fluid  saturated  porous  medium  [36]. 
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Stoll  and  Kan  [38]  derived  the  complex  reflection  coefficient  of  plane  acoustic 
waves  from  a  poro-elastic  sediment  half-space.  The  boundary  condition  model  is  based 
on  the  classical  work  of  Biot  that,  as  previously  stated,  predicts  three  different  body 
waves  in  the  sediment.  As  a  consequence,  when  plane  waves  in  water  are  incident  upon 
a  water-sediment  interface,  as  many  as  three  waves  can  be  generated  in  the  sediment,  de¬ 
pending  on  the  angle  of  incidence. 

Stoll  and  Kan  used  potentials  in  their  derivation  of  the  boundary  conditions. 
Johnson  and  Fiona  rewrote  Stoll  and  Kan’s  boundary  conditions  using  displacements  “as 
it  is  simpler  than  (although  equivalent  to)”  the  potential  derivation  [20].  Johnson  and 
Fiona’s  notation  for  the  boundary  conditions  produces  a  symmetric  determinant  when 
solving  for  the  eigenvalues  and  eigenvectors,  which  can  be  modeled  into  a  physically  re¬ 
alizable  system.  Also  the  density  terms  in  Johnson  and  Fiona’s  determinant  add  up  to  the 
density  of  the  water-saturated  sand.  Neither  of  these  are  true  for  Stoll  and  Kan’s  notation 
and  both  are  important  for  the  lumped  parameter  model  to  be  developed  here. 

1.4  OBSERVATIONS  OF  BIOT  WAVES 

Fiona  observed  a  second  bulk  compressional  wave  that  propagated  at  speeds  ap¬ 
proximately  25%  of  the  speed  of  the  normal  bulk  compressional  wave  in  a  fluid-saturated 
porous  medium  [29].  This  observation  was  made  using  an  ultrasonic  immersion  tech¬ 
nique  and  a  fluid-saturated  porous  medium  consisting  of  water  and  sintered  glass  spheres. 
The  excitation  of  the  slow  wave  in  the  solid  was  shown  to  be  consistent  with  the  princi¬ 
ples  of  mode  conversion  and  refraction  at  plane  liquid/solid  interfaces.  To  Fiona’s 
knowledge,  this  type  of  bulk  wave  had  not  been  previously  observed  at  ultrasonic  fre¬ 
quencies. 

Biot’s  theory  has  been  used  in  several  applications  of  acoustic  wave  propagation 
in  porous  media,  such  as  superfluid/superleak  systems  [18],  slow  waves  and  the  consoli¬ 
dation  transition  [19,  20],  gels  [21],  bones  [25],  porous  solids  [15,  23,  30],  marine  sedi¬ 
ments  [17,  35,  41,  43,  45],  pressure  diffusion  through  porous  media  [8],  and  air-saturated 
consolidated  porous  media  [16]. 
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Johnson  et  al.  [22]  observed  slow  compression  waves  in  consolidated  fluid- 
saturated  porous  media,  e.g.,  fused  glass  beads  (called  “Ridgefield  Sandstone”  even 
though  it  was  man-made)  and  ceramic  water  filters  (QF-20),  and  calculated  for  the  first 
time  all  of  the  input  parameters  necessary  for  a  complete  description  of  the  acoustic 
properties  of  all  the  modes  over  the  entire  frequency  spectrum. 

All  observed  slow  Biot  waves  have  been  in  consolidated  porous  media.  There 
have  been  no  uncontested  observations  of  a  slow  Biot  wave  in  an  unconsolidated,  e.g., 
sandy,  medium.  See  Controversy-2  that  follows. 

1.5  CONTROVERSY- 1  (Do  slow  Biot  waves  exist  in  sediments?) 

Even  though  Biot  theory  predicts  both  fast  and  slow  compressional  waves,  only 
the  fast  compressional  wave,  the  wave  of  the  first  kind,  and  the  shear  wave  have  been  ob¬ 
served  in  the  vast  majority  of  porous  systems.  Johnson  and  Fiona,  who  observed  slow 
Biot  waves  in  water-saturated  porous  media  where  the  solids  are  bonded  together,  have 
looked  for  slow  Biot  waves  in  unconsolidated  water  saturated  glass  beads  where  they 
conclude  “both  experimentally  and  theoretically  there  is  only  one  compressional  wave” 
[20]. 

Simpson  et  ah,  [34]  did  an  experiment  showing  the  wave  fields  in  sand  with 
smooth  and  rough  interfaces  at  all  angles  of  incidence.  They  concluded,  “This  data  does 
not  support  the  Biot  medium  model  predictions  of  a  second  slower  compressional  wave 
in  the  unconsolidated  water  saturated  porous  medium.” 

Seifert  et  ah,  [33]  did  an  experiment  with  sand  saturated  with  fluids  of  different 
viscosities  (water,  two  different  silicon  oils,  and  castor  oil)  to  measure  attenuation.  They 
stated,  “Biot’s  theory  shows  that  acoustic  waves  create  relative  motion  between  the  fluid 
and  the  solid  matrix  due  to  inertial  effects,  resulting  in  viscous  dissipation  of  acoustic  en¬ 
ergy.”  But  they  concluded  “The  attenuation  shows  no  correlation  with  the  viscosity  of 
the  different  pore  fluids  and  thus  theories  that  depend  upon  fluid  flow  cannot  explain 
these  data.”  They  conclude  Biot  theory  does  not  explain  attenuations  measured  in  fluid 
filled  sand. 
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On  the  pro  side,  Stoll  et  al.,  [38]  say  Biot  waves  exists  but,  “Waves  of  the  second 
kind  are  difficult  to  observe  in  sediments  where  water  is  the  saturant  because  of  their  high 
attenuation  and  the  way  in  which  energy  is  partitioned  at  the  interface.” 

1.6  CONTROVERSY-2  (Compressional  waves  at  1200  m/s  are  from  slow  Biot  or  scat¬ 
tering?) 

Recent  experimental  results  (Boyle  and  Chotiros  1992  [6],  Chotiros  1995  [9])  re¬ 
veal  acoustic  penetration  from  water  into  sandy  sediments  at  grazing  angles  below  the 
compressional  critical  angle  in  relation  to  the  mean  surface.  These  authors  interpret  their 
results  to  indicate  the  excitation  of  a  Biot  slow  wave  at  1200  m/s  in  the  sediment. 

But  numerous  authors  disagree.  Analytical  derivations  show  that  roughness  of 
the  water-sediment  interface  causes  propagation  of  acoustical  energy  from  water  into  the 
sediment  at  grazing  angles  below  the  compressional  critical  grazing  angle,  and  simula¬ 
tions  indicate  that  the  experimental  results  can  be  explained  in  terms  of  diffraction  of  an 
ordinary  longitudinal  wave.  By  modeling  sand  as  a  fluid,  and  including  a  small  amount 
of  roughness,  the  simulation  results  match  the  acoustical  penetration  experimental  results 
by  Boyle  and  Chotiros  (1992),  both  in  magnitude  and  in  arrival  time  [28,  40].  Mellema 
[27]  demonstrated  that  significant  energy  is  scattered  across  a  rough  fluid-solid  interface 
at  shallow  grazing  angles  and  that  the  received  intensity  can  be  accurately  modeled  using 
first  order  perturbation  theory.  For  the  roughened  case,  it  is  clear  that  the  later  time  arri¬ 
val  is  best  described  by  an  assumed  ray  path  in  the  water  to  a  point  above  the  buried  hy¬ 
drophone,  scattered  from  the  roughened  interface,  and  then  propagation  nearly  vertically 
to  the  buried  hydrophone  at  the  speed  of  the  ordinary  compressional  wave  (1680  m/s). 
This  can  have  the  appearance  of  a  wave  traveling  at  1200  m/s  [34]. 

Recent  field  measurements  [17,  26,  34]  have  strongly  supported  the  conclusion 
that  subcritical  penetration  in  sands  is  due  to  scattering,  not  Biot  slow  waves. 

To  match  the  1200  m/s  slow  compression  wave  using  Biot  theory,  Chotiros  [9] 
used  structural  and  material  moduli  of  sand  different  than  previous  authors.  Stoll  [39] 
and  other  authors  used  moduli  that  predict  a  slow  wave  having  a  speed  of  roughly  400 
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m/s.  Stoll  argues  Chotiros  is  not  justified  in  his  claim  that  a  Biot  wave  of  the  second 
kind,  with  a  wave  speed  of  1200  m/s,  had  been  detected  in  a  near-bottom  sand  deposit. 
This  opinion  is  based  on  the  fact  that  the  author  (Chotiros)  has  used  an  unacceptable 
value  for  the  bulk  modulus  of  the  skeletal  frame  and  a  questionable  value  for  the  bulk 
modulus  of  the  individual  grains.” 

Chotiros  [11]  replied  that  the  approach  advocated  by  Stoll  to  calculate  bulk 
moduli  works  well  for  a  porous  materials  with  cemented  frames,  such  as  fused  glass 
beads  and  sandstone,  but  it  does  not  work  for  sand.  Chotiros  also  contended  that  the 
measured  reflection  coefficient  for  normal  incidence  on  a  water-sand  interface  does  not 
match  elementary,  Rayleigh,  viscoelastic  reflection  theory  but  his  parameters  using  Biot 
theory  do  produce  a  match  with  the  measured  coefficient. 

1.7  PARAMETERS  AND  ASSUMPTIONS 

A  survey  of  the  literature  suggests  that  there  are  a  number  of  parameters  that  play 
a  principal  role  in  controlling  the  dynamic  response  of  saturated  sediments.  Of  these, 
Stoll  suggests  the  following  may  be  important  (not  necessarily  in  the  order  listed)  [37]: 

a.  dynamic  strain  amplitude, 

b.  porosity, 

c.  static  intergranular  stress, 

d.  gradation  of  sand  size  and  grain  shape, 

e.  material  properties  of  individual  grains, 

f  degree  and  kind  of  lithification  (bonding  of  particles), 
g.  structure  as  determined  by  the  mode  of  deposition. 

In  his  classic  development,  Biot  produced  constitutive  relationships  for  fluid- 
saturated  granular  media  and  followed  these  with  an  analysis  of  elastic  wave  propagation 
in  such  media  by  means  of  a  Eagrangian  formulation.  The  assumptions  underlying  this 
derivation  are  that  [1]: 

a.  The  medium  is  isotropic,  quasi-homogeneous  and  that  the  porosity  is  uniform 
throughout; 
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b.  The  pore  size  is  very  much  less  than  the  wavelengths  of  interest; 

c.  Scattering,  in  the  sense  of  diffraction  around  individual  grains  or  particles,  can  be 
ignored; 

d.  Pore  walls  are  impervious  and  the  pore  size  is  concentrated  around  an  average 
value; 

e.  Fluid  is  compressible; 

f  Fluid  may  flow  relative  to  the  solid. 

1.8  NEW  CONTRIBUTIONS  OF  THIS  DISSERTATION 

The  contributions  of  this  dissertation  are  both  theoretical  and  experimental.  On 
the  theoretical  side  the  primary  contributions  include:  (1)  Making  a  lumped  parameter 
model  from  the  continuous  parameter  Biot  theory,  fig.  2.1.  (2)  Deriving  the  fast  and  slow 
acoustic  impedances,  eqs.  (2-35)  and  (2-37).  (3)  Deriving  pressure  reflection  and  trans¬ 
mission  coefficients  from  impedances  and  effective  pressures,  tables  2.1  and  2.2.  (4)  De¬ 
riving  mode  shapes  from  effective  pressures  and  effective  densities,  eqs.  (4-27)  and  (4- 
28). 

On  the  experimental  side  the  primary  contributions  include:  (5)  A  reflection  ratio 
technique  in  which  reflection  amplitudes  are  divided  leaving  only  information  (pressure 
reflection  coefficient)  about  the  surface,  e.g.,  eq.  (4-1),  without  calibrating  the  transduc¬ 
ers  or  using  geometric  corrections.  (6)  Measuring  the  fast  and  slow  wave  acoustic  im¬ 
pedances,  fig.  4.1.  (7)  Measuring  the  fast  and  slow  effective  densities,  fig.  4.2.  (8) 
Bounding  the  slow  acoustic  wave  speed  in  an  unconsolidated  medium,  fig.  4.3.  (9) 
Measuring  the  fast  and  slow  wave  effective  pressures  from  an  open  pore  boundary  con¬ 
dition,  fig.  4.4.  (10)  Measuring  the  fast  and  slow  wave  mode  shapes,  fig.  4.6.  (11)  And 
last,  measuring  the  generalized  elastic  coefficients  (intermediate  moduli),  fig.  4.7. 
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Chapter  2: 

MATHEMATICAL  THEORY 

This  chapter  reviews  the  continuous  parameter  equations  of  motion  and  boundary 
conditions,  and  then  develops  a  lumped  parameter  model  from  which  the  fast  and  slow 
impedances,  moment  arms,  effective  densities,  and  effective  pressures  are  derived.  Then 
pressure  coefficients  are  derived  from  impedances  and  effective  pressures.  Only  nor¬ 
mally  incident  pressure  waves,  waves  propagating  perpendicular  to  the  sediment  bound¬ 
ary,  are  considered  in  the  following  continuous  parameter  review  and  lumped  parameter 
development.  This  consideration  simplifies  the  equations  by  removing  shear  waves  (and 
coupling  between  shear  and  pressure  waves),  and  is  sufficient  for  analysis  of  the  experi¬ 
mental  data  discussed  in  the  next  chapter. 


2.1  DIFFERENTIAL  EQUATIONS  OE  MOTION  AND  BOUNDARY  CONDITIONS 
Erom  Johnson  and  Fiona  [20],  eq.  (la)  and  (lb),  the  homogeneous  differential 
equations  of  motion  with  waves  propagating  in  the  z-direction  with  no  shear  displace¬ 
ments  are: 
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Eq.  (2-1)  is  the  equation  of  force  on  sand  particles  per  unit  area  (area  of  sand  plus  area  of 
fluid  in  sediment)  and  eq.  (2-2)  is  the  force  equation  on  fluid  particles  per  unit  area.  The 
unknowns,  Us  and  Uf,  are  displacements  (in  the  z-direction)  of  the  sand  and  fluid  particles. 
P,  Q,  and  R  are  the  generalized  elastic  coefficients  (intermediate  moduli). 
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Where  j8  is  the  porosity  (volume  of  fluid  per  unit  volume  of  sediment),  Kb  is  the  bulk 


modulus  of  the  sand’s  skeletal  frame,  Kg  is  the  bulk  modulus  of  the  sand,  K„  is  the  bulk 
modulus  of  the  water,  and  N  is  the  shear  modulus  of  the  skeletal  frame. 

The  density  terms  are: 
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Where  ps  is  the  density  of  the  sand  grains,  pf  is  density  of  water,  and  a  is  a  virtual  mass 


eonstant,  or  tortuosity,  expressing  the  apparent  inerease  in  the  mass  of  the  fluid.  As  a 
fluid  partiele  oseillates  in  the  z-direetion  it  must  translate  slightly  in  the  x-  and  y- 
directions  to  go  around  sand  particles  causing  it  to  appear  more  massive. 

The  viscous  term  is: 
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The  viscous  term,  F,  is  equivalent  to  bF,  eq.  (4.3)  in  [4],  where  b  is  defined  by  eq.  (6.8) 
in  [3].  In  eqs.  (2-9)  and  (2-10),  p  is  the  viscosity  of  water,  k  is  the  permeability  of  the 

medium,  Ji  and  Jo  are  Bessel  functions,  and  a  is  the  average  radius  of  the  pore  holes. 
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Solving  eqs.  (2-1)  and  (2-2)  for  plane  waves,  assume  the  solid  displacements  and 
fluid  displacements  are: 

=  y4exp(z(ft)t-fe))  (2-11) 

Uf  =  B  QXTp(i{cot  -  kz)) .  (2-12) 

Here,  k  is  wave  number,  t  is  time,  A  is  the  sand  displacement  amplitude  and  B  is  the  fluid 
displacement  amplitude.  Put  eqs.  (2-11)  and  (2-12)  into  (2-1)  and  (2-2)  and  express  the 
results  in  matrix  form: 
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The  square  matrix  on  the  right  hand  side  of  eq.  (2-14)  is  the  density  matrix  and  will  be 
referred  to  later.  Solve  for  the  eigenvalues  (wave  numbers  squared)  and  eigenvectors 
(mode  shapes).  This  will  give  displacements  of  the  sand  and  fluid  in  terms  of  the  wave 
numbers  and  the  mode  shapes. 

=  Ai  exp(z'(ft)t  -  kiz)^+  A2  exp(z(ft)t  -  k2z))  (2-15) 

Uj-  =  Bi  exp(z(ft)t  -  kiz)^+  B2  exp(z'(ft)t  -  k2z))  (2-16) 


Where  the  wave  fast  and  slow  numbers  are: 
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The  mode  shapes  are  (water  particle  displacement  divided  by  sand  particle  displace¬ 
ment): 
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The  boundary  conditions  for  the  water-sediment  (open  pore)  interface  are  taken 
from  Johnson  et  al.  [22]  and  specialized  to  normal  incidence  for  which  there  is  no  excita¬ 
tion  of  shear  waves.  Combine  equations  (2.17),  (2.18),  (2.19),  (2.20),  (2.21),  and  (2.25) 
from  ref.  [22]  and  get: 
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Eq.  (2-20)  says  the  force  in  the  sand  per  unit  area  of  interface  (note:  not  the  pressure  in 
the  sand  because  area  is  of  sand  and  pore  water)  plus  the  force  in  the  sediment  fluid  per 
unit  area  of  interface  is  equal  to  the  pressure  in  the  water  above  the  sediment.  Johnson’s 
W  is  the  same  as  this  dissertation’s  m,  +  Uy.  The  variables  m,  and  Uy  are  the  incident  and 
reflected  displacements,  assumed  to  be  of  the  form: 

Ui  =  exp(z(ft)t  -  k^z))  (2-21) 

M,  =  exp(z'(ft)t  +  k^z)) .  (2-22) 

Combine  equations  (2.18),  (2.20),  and  (2.23)  from  ref  [22]  and  get: 
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Eq.  (2-23)  says  the  force  in  the  sediment  fluid  at  the  interface  is  equal  to  the  force  in  the 
water  above  the  sediment  fluid  per  unit  area  of  interface.  Erom  equations  (2.22)  and 
(2.25)  from  ref  [22]: 

m^(1-;8)-i-m^j8  =  M;- -I-M, .  (2-24) 

Eq.  (2-24)  says  the  displacement  of  sand  times  the  area  of  sand  per  unit  area  of  interface 
plus  the  displacement  of  sediment  fluid  times  the  area  of  fluid  per  unit  area  of  interface  is 
equal  to  the  displacements  of  the  incident  plus  reflected  waves  in  the  water  per  unit  area. 
This  is  the  so-called  open  pore  boundary  condition.  Eq.  (2-24)  can  be  rewritten: 
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fast slow  (2-25) 

where 

u  fast  =  (4  (l  -  /^  )+  ^1  )exp(z(ry^  -k^z))  (2-26) 

Usiow  =  (^2(l-/^)+^2j8)exp(z(«^-^2^))-  (2-27) 

Eq.  (2-26)  is  the  average  of  the  sand  and  fluid  particle  displacements  during  the  passing 

of  a  fast  wave.  Both  previous  equations  will  be  used  later  to  derive  participation  factors. 
The  closed  pore  boundary  condition  for  displacements  will  be  used  later  and  is: 

u  =  Ui  +  Uj.  (2-27a) 

where  atz  =  0: 

u  =  u^  =  Uf  ={A^  +  A2  )exp(zft)t)  =  (fij  -I-  B2  )exp(zft)t)  (2-27b) 

This  section  (equations  of  motion  and  boundary  conditions)  has  given  the  wave 
numbers,  mode  shapes,  and  the  boundary  conditions  for  a  continuous,  2  degrees-of- 
freedom,  Biot  medium  having  a  planar  interface  and  normal  incidence.  The  next  section 
develops  the  lumped  parameter  model. 

2.2  LUMPED  PARAMETER  MODEL 

The  lumped  parameter  method  developed  here  was  strongly  influenced  by  the 
Dynamic  Design  Analysis  Method  (DDAM)  introduced  by  Belsheim  and  O’Hara  [2]. 
DDAM  is  used  to  analyze  non-contact,  bomb  shock  on  Naval  ships.  See  Appendix  E. 
DDAM  requires  a  physical  system  to  be  engineered  into  a  lumped  mass-spring  model. 
Then  the  equations  of  motion  are  written  in  matrix  form  from  the  mass-spring  model. 
The  natural  frequencies  and  mode  shapes  are  calculated  from  the  determinant  of  the  ei¬ 
genvalue  problem,  and  the  determinant  is  always  symmetrical.  Likewise,  even  though 
Biot’s  equations  of  motion  are  for  continuous  parameters,  the  determinant  for  the  wave 
numbers  is  symmetrical  using  Johnson  and  Plona’s  notation.  The  inspiration  for  the 
lumped  parameter  model  was  the  question  “Can  Biot’s  continuous  parameter,  symmetric 
determinate  be  reverse-engineered  (through  trial-and-error)  into  a  lumped  mass-spring 
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model?”  A  lumped  parameter  model  has  the  advantage  of  allowing  easy  visualization  of 
movements  and  drawing  free  body  diagrams  of  the  forces. 

The  continuous  medium  homogeneous  differential  equations  of  motion  and 
boundary  conditions  eqs.  (2-13),  (2-20),  (2-23),  and  (2-24)  can  be  described  by  a  lumped 
mass,  linear  spring,  concentrated  force,  viscous  damp-pot,  rigid-bar  mechanical  system, 
see  Fig.  2.1  and  Appendix  D. 


Figure  2.1 

Mass-Spring-Force  Mechanical  Model 
of  the  Biot  Medium 


In  going  from  the  continuous  model  to  the  lumped  model,  stiffness  (N/m)  and 
modulus  (Pa),  mass  (kg)  and  density  (kg/m^),  force  (N)  and  pressure  (Pa),  viscous 
damping  (Ns/m)  and  body  viscous  damping  (Ns/m"*),  difference  in  displacement  of  two 
points  (m)  and  strain  (m/m),  acoustic  impedance  (Pa-s/m)  and  mechanical  impedance 
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(Ns/m)  are  used  interchangeably.  Hopefully  this  does  not  cause  too  much  confusion. 
Point  H  in  Fig.  2.1  is  the  “location”  of  the  sand  (displacement  of  point  H  equals  the  sand 
particle  displacement)  and  point  B  is  the  “location”  of  the  fluid  in  the  sediment  (also  the 
displacement  of  point  B  equals  the  sediment  fluid  particle  displacement).  Pi  and  Pr  are 
the  incident  and  reflected  forces. 

To  show  that  the  boundary  conditions  of  the  mechanical  system  of  Fig.  2.1  are 
equivalent  to  Johnson  et  al.  [22]  boundary  conditions  draw  a  free  body  diagram  at  z  =  0, 
see  Fig.  2.2.  In  Fig.  2.2  the  positive  direction  is  down.  Eq.  (2-20)  is  the  sum  of  forces  in 
the  vertical  direction.  Eq.  (2-23)  is  the  sum  of  moments  about  point  H.  The  incident  and 
reflected  displacements,  m  and  Ur,  are  related  to  the  solid  and  fluid  displacements,  Us  and 
Uf,  by  eq.  (2-24).  Eig.  2.3  is  equivalent  to  Eig.  2.2  but  in  Rayleigh  notation,  i.e.,  e'®^  sup¬ 
pressed  (ref  [31],  section  270,  eq.  (11)). 


z  =  0 


i(P  -  Q)(kiA]  +  k2A2) 


i2Q(kiAi  +  k2A2  ... 
+  kiB]  +  k2B2) 


- > 

i(R  -  Q)(k,Bi  +  k2B2) 


Figure  2.3 

Rayleigh  Notation  of  Figure  2.2 


Figure  2.4 

Free-Body  Diagram  of  Fast  Wave  Mode  Shape 
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A  force  can  be  positioned  that  will  only  excite  the  fast  wave.  As  will  be  seen,  this 
is  a  consequence  of  mode  orthogonality.  See  Fig.  2.4.  From  Fig.  2.4,  sum  forces  in  the 
vertical  direction: 

Fj  =  iiPk^A^  +  Qk^A^  +  QkiBi  +  Rk^B^).  (2-28) 

The  sum  of  the  moments  around  point  A  in  Fig.  2.4  is: 

<5iFi  =  iiQkiA^  +  RkiBi ) .  (2-29) 


From  equations  (2-28)  and  (2-29)  the  moment  arm  for  a  force  to  excite  only  mode  1  is: 


QAi  -I-  RB^ 

PAi+  QAi  +  QB^  +  RBi 


(2-30) 


Another  way  to  find  (5i  is  by  knowing  that  a  foree  applied  at  82,  the  moment  arm 


for  mode  2  (slow  wave),  will  cause  no  work  at  8^,  i.e.,  the  translation  at  8^  in  the  z- 
direction  is  zero.  See  Fig.  2.5. 

A. 

(2-31) 


(5i  = 


^2  ~  Pi 


Fi 


-A 


^  ^A 


B 


B 


Fy 


Figure  2.5 

Location  of  F\  Using  Orthogonality 


These  two  relationships  for  8^,  equations  (2-30)  and  (2-31),  are  equivalent  to  the 


orthogonality  relationship: 
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jA 

\b. 


P  Q 
Q  R 


\A 

\B, 


=  0. 


From  Fig.  2.4  the  velocity  at  point  5^  is: 

vj  =z(l-<5i)^^+z(5i  ^;^  =  z[(l-<5i)m4  +5y(oB{\. 
at  at 

The  impedance  of  first  mode  (fast  wave)  is: 

Z.=^ 

^1 


21  = 


ki  [(P  +  Q)A,+(R  +  Q)bJ 
(O  +2QA^B^+RB]^ 


(2-32) 


(2-33) 

(2-34) 

(2-35) 


Similarly: 


QA2  +  RB  2 


A, 


PA2  +  QA2  +  QB2  +  RB2  Ai  -  Bi 

^  _k2  [iP+Q)A2  +  iR  +  Q)B2l 

^  0)  PA2^  +  2QA2B2  +  RB2^ 


Because: 


k 


f  Z.A2 


K, 


(2-36) 

(2-37) 

(2-38) 

(2-39) 


the  modal  stiffnesses  (moduli)  and  masses  (densities)  of  the  first  and  second  modes  are: 

[{P  +  Q)A^+{R  +  Q)bJ 


Ki  = 


K2  = 


Pi 


PA^^  +2QA^B^+  RB^^ 

[{P  +  Q)A2+iR  +  Q)B2l 

PA2^  +2QA2B  2+ RB  2^ 

[iP  +  Q)A,+{R  +  Q)Bj 
PA^^  +2QA^B^+RBi^ 


(2-40) 

(2-41) 

(2-42) 
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P2  - 


^k2f[iP  +  Q)A2+iR  +  Q)B2f 

(0  )  PA2^+2QA2B2+RB2^ 


(2-43) 


The  following  two  equations  are  taken  from  the  matrix  equation  of  (2-14): 


(PA  +  QB)-- 


(QA  +  RB)-- 


f 

( 

Pi  1  “ ' 

A  + 

P12  +  '~ 

-V 

CO  J 

1 

f] 

( 

- 

P12  + 

— 

A  + 

P22  — 

B 

x 

CO  ^ 

1 

_ 

(2-44) 

(2-45) 


Add  equations  (2-44)  and  (2-45)  and  put  in  the  numerators  of  (2-42)  and  (2-43).  Also 
multiply  equation  (2-44)  by  A  and  (2-45)  by  B,  add,  and  put  in  the  denominators  of  (2- 
42)  and  (2-43).  The  results  are  the  modal  masses  or  the  effective  masses  in  terms  of  den¬ 
sity  (and  viscous)  parameters: 

[(Pi  1  +  Pi2  )Ai  +  (Pi2  +  P22  )^i  ? 


Pi  = 


F 


Pi  1^1  +2pi2^1^1 +P22^1  (4l“^l) 


CO 


P2  - 


[(Pl  1  +  P12  )^2  +  (Pl2  +  P22  )^2  P 


P\\A2^  +2pi2^2^2  +  P22^2^  — (^2  “  ^2  P 

CO 


(2-46) 


(2-47) 


These  modal  masses,  eqs.  (2-46)  and  (2-47),  agree  with  the  Dynamic  Design 
Analysis  Method  (DDAM),  see  Appendix  E.  Also,  Clough  and  Penzien  [14]  p.  559, 
show  that  the  sum  of  the  modal  masses  is  equal  to  the  sum  of  all  the  terms  in  the  mass 
matrix,  eq.  (2-14),  i.e.,  pi  -I-  p2=  pii  +  P22+  2pi2=  (l-j8)Pi  +  Ppf=  p.  Where  pis  the  total 

density  of  the  water  filled  sand.  With  Johnson  and  Fiona’s  form  of  the  Biot  equations, 
eqs.  (2-1)  and  (2-2),  the  modal  masses  add  up  to  the  total  mass,  i.e.,  pi  +  p2=  p. 


The  terms  that  link  particle  displacements  to  effective  (or  modal)  displacements 
are  called  “participation  factors”  in  this  dissertation.  For  a  single  degree-of-freedom 
system,  of  course,  particle  displacements  and  effective  displacements  are  equal.  Trans¬ 
forming  average  particle  displacement  of  the  sand  and  fluid  during  the  passing  of  a  fast 
wave,  eq.  (2-26),  into  a  lumped  parameter  model.  Fig.  2.6,  shows  the  average  particle 
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displacement  is  located  at  a  distance  from  point  A.  In  Fig.  2.6  points  1  and  2  are  the 

locations  of  the  first  and  second  modes  and  points  A  and  B  are  the  locations  of  the  sand 
and  fluid  particles.  The  symbols  ui  and  U2  are  the  effeetive  or  modal  displacements.  The 
fast  wave  participation  factor  is: 


p  _  ^fast 
M  - - 

Ui 

82 

And  likewise  the  slow  wave  participation  factor  is: 

p  _  ^slow 

^2 - 

U2 

_  P-s, 

82  -81 


(2-48) 


(2-49) 
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Fig.  2.7  is  a  mechanical  system  equivalent  to  Fig.  2.1  but  in  terms  of  modal  pa¬ 
rameters.  The  acoustic  impedances  for  Biot  waves  are  fully  described  by  the  impedances 
and  the  impedance  moment  arms,  which  are  modal  parameters  independent  of  the  load¬ 
ing  of  the  mechanical  system,  i.e.,  they  are  homogeneous  parameters  (in  a  differential 
equation  sense).  The  terms  that  link  loading  (boundary  conditions)  and  modal  parame¬ 
ters  are  the  modal  forces  or  effective  forces. 


•  • 

•  • 

*  Figure  2.7  * 

Mechanical  System  Equivalent  to  Fig.  2.1 
But  In  Terms  of  Modal  Parameters 

In  the  open  pore  boundary  mentioned  earlier,  the  water  above  the  sediment  and 
the  fluid  in  the  sediment  are  in  direct  contact.  This  is  the  most  commonly  used  boundary 
condition  for  sediments.  The  terms  that  link  external  forces  to  forces  on  the  first  and 
second  mode  are  called  “effective  forces.”  The  incident  and  reflected  forces  act  at  a  dis¬ 
tance  (5  from  point  A  for  an  open  pore  interface,  see  Fig.  2.7.  Let  the  incident  and  re- 
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fleeted  forces  equal  one  then  from  Fig.  2.8,  P\  is  the  “effective  force”  on  mode  1  and  P2 
is  the  “effective  force”  on  mode  2. 


«2-«l 

Si-S, 


(2-50) 

(2-51) 


Note  that,  numerically,  open  pore  effective  forces  are  equal  to  participation  fac¬ 
tors  even  though  their  meanings  are  different. 

For  a  closed  pore  boundary  the  sediment  is  in  contact  with  a  non-porous  material, 
like  acrylic.  In  a  closed  pore  boundary  the  sand  and  fluid  particles  and  the  fast  and  slow 
modal  displacements  move  with  the  same  amplitude  at  the  interface,  and  the  modal 
points  1  and  2  move  equally,  see  eq.  (2-27b).  The  forces  needed  to  move  these  points 
equally  are  the  closed  pore  effective  forces. 


(,)  ^  ^1 

1  Z1+Z2 


(2-52) 


)(c) 


Zi  +Z2 


(2-53) 


The  next  section  calculates  the  pressure  reflection  and  transmission  coefficients 
using  impedances  and  effective  forces. 
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2.3  REFLECTION  AND  TRANSMISSION  PRESSURE  COEFFICIENTS  FOR  A 
SOUND  WAVE  IN  WATER  INCIDENT  UPON  A  SAND  INTERFACE  (OPEN  PORE) 
This  section  finds  the  reflection  and  transmission  coefficients  for  various  bound¬ 
ary  conditions  of  a  medium  that  supports  Biot  waves.  If  the  distance  from  the  first  to 
second  modes,  points  I  and  2,  are  normalized  to  one,  the  normalized  distances  from  an 
applied  unit  force  (the  force  is  applied  at  a  position  to  make  the  forces  on  Z\  and  Z2  equal 
to  P\  and  Pi)  to  the  locations  of  the  modes  are  numerically  equal  to  the  normalized  effec¬ 
tive  forces,  P\  and  P2,  see  Fig.  2.9.  These  normalized  forces  and  normalized  distances 
will  be  assigned  the  same  symbols,  P\  and  P2.  The  impedance  of  the  Biot  medium  is 
found  by  applying  a  force  of  1  at  a  distance  P2  from  the  left  in  Fig.  2.9.  The  springs  in 
the  figure  are  labeled  as  impedances,  Z  (units  of  Ns/m  in  lumped  system).  Technically, 
the  spring  constants  should  be  Zco  (N/m),  but  the  frequency  terms  are  dropped.  Veloci¬ 
ties  of  the  first  and  second  modes  are: 


V,  = 


Vt  =■ 


(2-54) 


The  velocity  at  the  point  of  unit  force  is: 


v  =  viE)  +V2P2  = 


Z^Pi'+ZiP;- 

Z,Z, 


The  impedance  “seen”  by  an  incident  wave  from  water  to  a  Biot  medium  is: 


z  =  l 

V 


Z  = 


Z1Z2 


Z2V+Z1E2' 


(2-55) 


(2-56) 


The  pressure  reflection  coefficient  of  a  sound  wave  off  a  Biot  medium  (water  to 
sand)  is  (where  Z^  is  the  impedance  of  water): 

Z  -Z„, 


Z  +  Z., 


Z1Z2  -I-  Z^Zji2  + 


(2-57) 
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The  pressure  transmission  eoeffieient  for  the  fast  and  slow  waves  eombined  is: 


2Z 

^ 


T'.s  = - ^  •  (2-58) 

ZiZ2+Z^Zii2 

From  the  free-body  diagram  of  Fig.  2.10  the  fast  wave  transmission  eoeffieient  is: 
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^  w\ 


T  =  PT 

w\  ■'  H  wi" 


2P1Z1Z2 


Z1Z2  +  Z^Zj/2  + 


And  the  slow  wave  transmission  coeffieient  is: 

T  =  P-,T 

^  w2  ^  2-^  ws 


Twi  - 


2P2Z1Z2 


Z1Z2+  Z^ZiP2  +Z^Z2-Pi 


(2-59) 


(2-60) 


2.4  TRANSMISSION  AND  REFLECTION  COEFFICIENTS  OF  A  FAST  WAVE 
FROM  SAND  TO  WATER  (OPEN  PORE) 

The  impedanee  diagram  for  a  fast  wave  in  sand  ineident  on  a  water  interfaee  is 
shown  in  Fig.  2.11.  The  water  partiele  velocity  at  a  point  P2  from  the  left  from  unit  force 
at  the  left  is: 


= 


1  1 
Pi 


The  slow  wave  particle  velocity  is: 


^J_ 

Pi  Z2 


Therefore  the  fast  wave  particle  velocity  from  a  unit  force  at  the  fast  wave  is: 


1  Pi 

1  w  2 


^  1 

\  ^2" 

—  -1- 

[Plj 

7 

The  impedance  of  a  fast  wave  incident  on  a  sand-water  interface  is: 


(2-61) 


(2-62) 


(2-63) 


Z  -  1  - 
^Iw  -  -  ■ 


PlZ^Z2 


N  Z2  +P2  Z^ 


(2-64) 


The  transmitted  pressure  coefficient  from  a  fast  wave  into  the  water  and  into  the  slow 


P1-W&2 


Ziw+Zi 


wave  is: 


26 


'  1-W&2 


2^1 


Z1Z2  +  Pi  Z^Z2  +  P2 


(2-65) 


From  Fig.  2.12  the  transmitted  wave  into  the  water  from  an  ineident  fast  wave  is: 

Tlw  =~Ti-w&.2 

Ti^  = - ^ -  (2-66) 

-^1-^2  +^1  +^2 

The  refleeted  pressure  eoefficient  (refleeted  fast  wave  from  an  ineident  fast  wave,  open 
pore,  water  interfaee)  is: 

p(w)  _  -^Iw  ~  -^1 
11  7  j.  7 


P\^w^2  ^^\^2  ~P2^w^\ 
Pi  ^w^2  +  2^1 2^2  P2  '^w^l 


(2-67) 
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T  —  —T 

^\w  ^  1-W&2 


(Transmitted) 


Pi 


1  (Incident,  Fast) 


water 


-02 


sedi- 

ment 


A 

(Reflected  Fast) 


‘  1-W&2 


(Reflected  Slow) 


Figure  2.12 

Free-Body  Diagram  of  Incident  and  Reflected 
Forces  on  Sand-Water  Interface 


The  reflected  slow  wave  pressure  coefficient  is: 


dCw)  _  ^2  rp 
^12 


n(w) 

^12 


^P\P2^w^2 


Z1Z2  -I-  Pi  Z^Z2  -I-  P2  Z^Zi 


(2-68) 


2.5  TRANSMISSION  AND  REFLECTION  COEEEICIENTS  OE  A  SLOW  WAVE 
EROM  SAND  TO  WATER  (OPEN  PORE) 

The  impedance  of  a  slow  wave  incident  on  a  sand-water  interface  is: 


Z 


2w 


P2^Z^Zi 

Z\  +P^Z^ 


(2-69) 


The  transmitted  pressure  coefficient  from  a  slow  wave  into  the  water  and  into  the  fast 


wave  is: 
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2Z 


'  2-W&1 


2w 


^2w  +2^^ 


'  2-W&1 


^P2Z,Z, 


Z,Zo  +P?Z^.,Z.  +  A^Z,„Z, 


"1^2  ^w^2  "^^2 

The  transmitted  wave  into  the  water  from  an  ineident  slow  wave  is: 

^2w  =  ~^2-w&l 
^2 


(2-70) 


^2w  - 


Z1Z2  +  Pi  Z^Z2  +  P2 


(2-71) 


The  refleeted  pressure  eoeffieient  is  (refleeted  slow  wave  from  an  ineident  slow  wave, 
open  pore,  water  interfaee): 


^2w  ^2 

-^2w  +  -^2 


p(w) 

^22 


^\^2  P\^w^2 
P\  Z,„Z7  -I-Z1Z7  +  ^2 


1  ^w^2  ■^'^1^2 

The  refleeted  fast  wave  pressure  eoeffieient: 

p('^)  _  T' 

^21  -^^2-w&l 

^2 


n(w) 

^21 


Z1Z2  +  Pi  Z^Z2  +  P2  Z^Zi 


(2-72) 


(2-73) 


2.6  REFLECTION  AND  TRANSMISSION  COEFFICIENTS  OF  A  FAST  WAVE 
FROM  SAND  TO  NON-POROUS  (CLOSED  PORE)  INTERFACE 

Figure  2.13  is  the  impedanee  diagram  for  a  fast  wave  ineident  on  a  non-porous 
interfaee  that  has  a  elosed  pore  surfaee  (e.g.,  sand  on  aerylie,  sand  on  Ethafoam,  ete.). 
The  sand  partieles,  water  partieles,  first  mode,  and  seeond  mode  all  have  the  same  dis- 
plaeement  on  a  closed  pore  surface;  therefore  the  fast  and  slow  wave  particle  velocities, 
V]  and  V2,  are  the  same  at  the  interface.  Impedance  seen  by  fast  wave  (Zc  is  impedance  of 
a  closed  pore,  non-porous  medium): 
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The  slow,  pressure  refleetion  eoeffieient  from  a  fast  ineident  wave  off  a  elosed  pore  sur- 
faee  is: 


30 


= 

^12 


-2Z. 


Zt  +  Z  „  +  Z 1 


(2-77) 


2.7  REFLECTION  AND  TRANSMISSION  COEFFICIENTS  OF  A  SLOW  WAVE 
FROM  SAND  TO  NON-POROUS  (CLOSED  PORE)  INTERFACE 
The  impedance  seen  by  slow  wave  is: 

Z  =  Zj  +  Z^  .  (2-78) 

The  slow,  pressure  reflection  coefficient  from  a  slow  incident  wave  off  a  closed  pore  sur¬ 
face  is: 

Z-Z. 


R 


(c) 
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Z  +  Zj 


(c)  _ 


z,  +z. 


R^  ’  = 

Z2+Z, +Zi 

The  transmitted  coefficient  from  a  slow  incident  wave  on  a  closed  pore  surface  is: 

2Z. 


(2-79) 


'2c 


Z  +  Z2 

2Z. 


'2c 


^2  +-^c  +-^1 


(2-80) 


The  fast,  pressure  reflection  coefficient  from  a  slow  incident  wave  off  a  closed  pore  sur¬ 
face  is: 


R^‘^^  = 
^21 


-2Z 


Z  +  Z 


r{c) 

^21 


■2Z, 


Z2  -l-Z^  -l-Zi 


(2-81) 
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2.8  REFLECTION  AND  TRANSMISSION  COEFFICIENTS  OF  AN  INCIDENT 
WAVE  FROM  A  NON-POROUS  MATERIAL  TO  SAND 

Figure  2.14  is  the  impedanee  diagram  for  a  fast  wave  ineident  on  sand  from  a 
elosed  pore  surfaee.  The  impedanee  seen  by  the  ineident  wave  is: 

Z  =  Zj  +  Z2 .  (2-82) 

The  transmission  of  the  fast  wave  is: 


■^1  +-^2 
-^1 


(2-84) 


(2-85) 


1 


Figure  2.14 
Ineident  Wave 

_ from _ 

Non-Biot  Material 
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2.9  SUMMARY  OF  PRESSURE  COEEEICIENTS 

The  Tables  2.1  and  2.2  summarize  the  pressure  eoefficients.  The  energy  balances 
for  these  coefficients,  i.e.,  pressure  times  velocity  of  the  incident  wave  equals  pressure 

times  velocity  of  the  response  waves,  e.g.,  —  =  H — iti  Also  the  pres- 

Zi  Zj  Zj 

sures  on  each  side  of  the  interface  are  equal,  i.e.,  one  plus  reflected  equals  transmitted, 
e.g.,  1  +  =  Ti^  (given P\+  ?!=  1).  Appendix  B  numerically  checks  parts  of 


Table  2.1.  The  closed  pore  table  (Table  2.2)  can  also  be  calculated  from  the  open  pore 
table  by  using  the  effective  pressures  for  closed  pore,  eqs.  (2-52)  and  (2-53). 

Pressure  coefficients  are  used  in  the  next  chapters  where  different  combinations 
of  acoustic  measurements  are  taken  to  extract  impedances,  effective  densities,  wave 
speeds,  effective  pressures,  mode  shapes,  and  intermediate  moduli.  Only  four  of  the  18 
pressure  coefficients  {Rcc,  T\c,  I’d,  and  Rn)  derived  are  used  in  this  experiment.  The 
other  14  coefficients  are  included  only  for  completeness. 


Table  2.1  -  Pressure  Coefficients,  Open  Pore,  Sand-Water  Interface 
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Chapter  3: 

ACOUSTIC  MEASUREMENTS 

The  Applied  Physies  Laboratory  has  an  acoustic  water  tank  system  suitable  for 
measuring  pressure  amplitudes  of  waves  transmitted  and  reflected  through  various  me¬ 
dia.  In  this  chapter,  specifications  are  reviewed,  reflection  equations  are  derived,  and 
frequency  measurements  are  discussed. 

3.1  SPECIFICATIONS 

The  water  tank  is  1.6  m  deep,  2.3  m  diameter,  and  holds  5000  liters  of  water.  The 
tank  has  a  lid  that  can  be  sealed  and  a  vacuum  pump  can  be  used  to  de-gas  the  contents  in 
the  tank.  The  tank  has  a  piping  and  filter  system  that  can  recirculate  and  purify  the  wa¬ 
ter.  A  resistivity  probe  monitors  purity  and  indicates  when  the  filters  need  to  be  changed. 
With  the  de-gassing  and  purifying  systems  the  water  is  “pure”  and  the  sand  is  bubble 
free.  The  tank  has  a  positioning  system  for  the  transmitter  and  receiver  with  an  accuracy 
of  about  25  |im  and  a  computer  based  measuring  system  including  a  data  acquisition  sys¬ 
tem  and  waveform  generator.  The  data  acquisition  employs  a  12  bit  analog- to-digital 
converter  and  samples  at  10  MHz.  An  overhead  electric  crane  handles  all  heavy  items 
going  in  and  coming  out  of  the  tank. 

The  pulses  were  shaped  so  that  the  residual  ringing  of  the  transducer  was  actively 
canceled  (see  Appendix  G).  Return  pulses  were  windowed  with  a  Gaussian  window  30 
|is  wide.  Measurements  of  each  surface  were  taken  at  100  different  locations  with  100 
samples  at  each  location.  The  raw  data  shown  in  Figs.  3.2,  3.4,  3.6,  and  3.8  are  the  co¬ 
herent  time  averages  of  100  samples  at  100  locations  (10,000  total  samples). 

The  sand  used  was  Ottawa  sand,  with  an  average  size  of  279  |im  and  a  standard 
deviation  of  about  100  |im.  Appendix  H  has  a  more  detailed  analysis  of  sand  size  and 
composition. 
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3.2  REFLECTION  EQUATIONS 

A  reflection  ratio  technique  was  used  to  extract  desired  information.  The  strategy 
behind  the  reflection  ratio  method  is  to  form  a  ratio  that  removes  most  unknowns  from 
measured  reflections.  Measurements  were  taken  so  that  two  different  acoustic  pulses  had 
exactly  the  same  propagation  path  but  different  reflection  surfaces  so  that  dividing  the 
two  acoustic  pulses  factored  out  all  unknowns  in  the  transmitter’s  signal,  the  acoustic 
path,  and  the  hydrophone’s  reception,  leaving  only  information  about  the  surfaces. 
Jumping  ahead  a  little,  if  the  reflection  from  the  second  surface  of  acrylic-sand,  r2s,  in 
Fig.  3.1  is  divided  by  the  reflection  of  the  second  surface  of  acrylic-air,  r2air,  in  Fig.  3.3 
what  remains  is  only  information  about  the  second  surface.  Mathematically  eq.  (3-3)  is 
divided  by  eq.  (3-2)  and  what  is  left  is  eq.  (4-1),  which  is  greatly  simplified  and  has  only 
two  unknowns.  These  reflection  ratios  were  then  used  to  find  selected  acoustic  parame¬ 
ters  in  Chapter  4.  To  this  end,  sand  was  held  in  a  tray  where  the  transmitter  and  receiver 
could  be  positioned  below  the  tray.  Then  the  top  surface,  of  the  sand  or  the  tray,  could 
be  changed  without  disturbing  the  path.  Plans  of  the  tray  and  lid  are  shown  in  Appendix 
F. 

The  analysis  in  Chapter  4  shows  there  are  14  unknowns  with  4  equations  from 
reflection  ratios  and  10  other  supporting  equations.  The  4  equations  from  reflection  ra¬ 
tios  were  taken  from  6  different  interfaces,  which  are:  (I)  water-acrylic,  (2)  acrylic-air, 
(3)  acrylic-sand,  (4)  sand-air,  (5)  sand-Ethafoam,  and  (6)  sand-water. 

The  complex  reflection  amplitude  from  the  first  surface  of  acrylic  (water-acrylic) 
is  (e.g.,Fig.  3.1): 

Rxa  =  (3-1) 

R\a  is  the  Fourier  Transform  of  (Fig.  3.2),  At  is  the  pressure  amplitude  a  unit  distance 
from  the  acoustic  center,  whose  distance  from  the  interface  is  dia/2.  Za  is  the  impedance 
of  the  4-inch  acrylic,  is  the  impedance  of  water,  k„  is  the  wave  number  of  water,  and 
is  the  distance  the  sound  travels  through  water  to  the  acrylic  surface  and  back.  Due 
to  the  mechanics  of  the  transducers  d\a  may  not  equal  dJ^^K  Thorsos  et  al  [40]  shows  eq. 
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(3-1)  is  a  good  approximation  when  is  much  larger  than  the  wavelength 

approximately  0.127  m,  wavelength  approximately  0.006m). 

Assuming  the  impedance  of  air  is  very  small  compared  with  acrylic,  the  reflection 
from  the  second  surface  with  air  (acrylic-air)  is  (Fig.  3.3): 
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(3-2) 
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Where  d2  is  a  transmission  loss  factor  and  can  be  considered  as  twice  the  distance  from 
the  “virtual”  acoustic  center  to  the  second  acrylic  surface.  This  is  not  equal  to  the  dis¬ 
tance  from  the  transducer  to  the  second  acrylic  surface  because  the  propagation  through 
two  media  forms  a  different  focusing  point.  The  wave  number  of  the  acrylic  is  ka,  and 
the  thickness  of  the  acrylic  is  ta. 

The  reflection  from  the  second  surface  with  sand  (acrylic-sand)  is  (Fig.  3.5): 
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Where  {RcA  is  the  reflection  of  a  wave  from  a  closed  pore  material  back  to  the  closed 
pore  material  from  a  closed  pore-sand  interface  where  acrylic  is  the  closed  pore  material. 
Rcc  is  taken  from  Table  2.2. 

The  reflection  of  the  third  surface  with  water  (sand-water)  is  (Fig.  3.1): 
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Where  dj,  is  a  transmission  loss  factor,  (Tc\)a  is  the  transmission  coefficient  of  a  wave 
from  a  closed  pore  material  to  a  fast  wave  in  sand  where  the  closed  pore  material  is 
acrylic,  (T \c)a  is  the  transmission  coefficient  of  a  fast  wave  from  sand  to  a  closed  pore 
material  where  the  closed  pore  material  is  acrylic,  RiA^  is  the  reflection  of  a  fast  wave 
back  to  a  fast  wave  from  a  sand- water  interface.  I’d  and  Tic  are  taken  from  Table  2.2. 
7?  11*^^  is  taken  from  Table  2.1. 

The  third  surface  with  Ethafoam  (sand-Ethafoam)  is  (Eig.  3.5): 
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Where  (Ri/‘^^)e  is  the  refleetion  of  a  fast  wave  to  a  fast  wave  from  a  sand-elosed  pore  in- 
terfaee  where  Ethafoam  is  the  elosed  pore  material.  is  taken  from  Table  2.2. 

The  third  surfaee  with  air  (sand-air)  is  (Fig.  3.7): 
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Where  4  is  the  thiekness  of  the  sand. 

Refleetion  ratios  are  obtained  with  the  lid  (2  ineh  aerylie  and  Ethafoam)  on  and 
off.  The  weight  of  the  lid  underwater  with  Ethafoam  attaehed  is  2.641  kg.  The  eontaet 
area  of  the  lid  is  0.4560  m  .  This  gives  an  overburden  pressure  equivalent  to  a  depth  of  5 
mm  of  water  saturated,  buoyant  sand.  The  resulting  small  stress  should  not  influenee  the 
refleetion  ratios. 


3.3  FREQUENCY  MEASUREMENTS 

In  the  following  ehapters  the  frequeney  domain  refleetions,  R,  are  the  Fourier 
transforms  of  the  time  domain  refleetions,  r.  The  refleetion  data  for  eaeh  loeation  were 
averaged  eoherently  in  the  frequeney  domain  of  100  transmissions. 
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Figure  3.2 

Sand  -  Water  Interface  {r^„) 
Acrylic  -  Sand  Interface  (ri^) 
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Figure  3.3 

Schematic  for  Acrylic-Air  Reflection  (r2a,>) 
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Figure  3.4 

4”  Acrylic  -  Air  Interface  (r2a;>) 
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Figure  3.5  ^ 

Schematic  for  Sand-Ethafoam  Reflection  {r^E) 
and  Water-Acrylic  Reflection  (ria) 
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Figure  3.6 

Sand  -  Ethafoam  Interface  (rs^) 
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Figure  3.7 

Schematic  for  Sand-Air  Reflection  {r^ait) 
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Figure  3.8 

Sand  -  Air  Interface  {r^air) 
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Chapter  4: 

ANALYSIS 

This  chapter  uses  the  theory  in  Ch.  2  and  the  measurements  in  Ch.  3  and  solves 
for  various  physical  parameters  and  acoustical  properties,  which  include:  fast  and  slow 
impedances,  fast  and  slow  wave  speeds,  fast  and  slow  effective  densities,  fast  and  slow 
effective  pressures,  fast  and  slow  mode  shapes,  the  three  intermediate  moduli,  and  a  term 
that  will  be  called  the  mass  factor,  which  is  a  combination  of  viscosity  and  tortuosity 
terms.  But  the  measurements  for  three  other  parameters,  the  sand  frame  modulus,  sand 
particle  modulus,  and  shear  modulus,  were  not  made  successfully  and  the  reasons  are 
discussed.  Other  parameters  that  were  measured  for  this  chapter’s  analysis  include:  the 
density  of  the  water  saturated  sand,  the  density  of  the  sand  particles,  the  impedance  of  the 
acrylic,  the  impedance  of  the  Ethafoam,  the  thickness  of  the  sand,  the  porosity  of  the 
sand,  and  the  volume  of  the  water  saturated  sand. 

4.1  FAST  AND  SLOW  WAVE  IMPEDANCES 

Reflection  ratios  are  used  to  calculate  the  fast  and  slow  wave  impedances.  The 
use  of  reflection  ratios  uses  relative  measurements  and  removes  the  need  for  measuring 
absolute  phases  and  amplitudes  (calibrating  the  transducers  and  compensating  for  propa¬ 
gation  effects).  From  the  reflections,  the  fast  and  slow  wave  impedances  are  calculated 
by  the  following  formulas.  The  pressure  reflection  coefficient  ratio  of  the  second  surface 
is  (combine  eqs.  (3-3),  (3-1),  and  (3-2)): 

i^2s  )// 
i^la  )// 

(^2air  )/ 

{Rla  )l 

The  subscripts  /  and  II  represent  acoustic  pulses  in  which  the  transducer  generated  ampli¬ 
tude,  At,  may  be  slightly  different.  R2air  and  R2s  are  divided  by  R\a  to  normalize  the  am¬ 
plitudes.  The  acoustic  centers  and  distance  traveled  through  water  were  held  the  same. 
Unfortunately  phase  measurements  were  not  accurate.  That  is  why  the  absolutes  of  the 
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reflections  and  the  negative  in  front  of  the  right  hand  side  are  used  in  eq.  (4-1).  It  is  felt 
the  phases  were  not  accurate  because  of  possible  disturbances  of  the  top  sand  surface 
(Ethafoam  on  and  off).  The  wavelength  of  sound  in  sand  was  6.6  mm  so  a  0.55  mm 
change  in  sand  thickness  (travel  distance  of  1.1  mm)  would  make  a  60°  difference  in 
phase.  The  error  from  taking  absolute  reflections  is  estimated  in  Appendix  C.  Also  the 
pressure  reflection  coefficient  of  the  third  surface  is  (combine  eqs.  (3-5),  (3-1),  and  (3- 
6)): 
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where  Zg  and  Za  are  measured  in  Appendix  A.  From  eqs.  (4-1)  and  (4-2),  the  imped¬ 
ances  of  the  fast  and  slow  waves  are: 

^ E  ^3E _3air  )(^  ^2s _2air  ~  ^2s _2air  )(^  ^3E _3air  ) 

2(1 +  ^25  2air  ) 


2l  = 
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(4-4) 


The  fast  and  slow  impedances  were  measured  in  100  different  locations.  The  av¬ 
erage  fast  and  slow  impedances  and  their  standard  deviations  are  shown  in  Fig.  4.F 

The  bands  or  tolerances  on  the  following  plots  are  all  sample  standard  deviations 
and  not  expected  error.  At  band  center  (250  kHz)  the  measurement  errors  (Appendix  C) 
are  about  one-third  as  large  as  the  observed  fluctuations  in  measurements;  therefore  this 
fluctuation  must  be  ascribed  to  spatial  variability  of  the  sand  sample.  At  the  band  edges 
(150  kHz  and  350  kHz,  the  pulse  has  a  bandwidth  of  200  kHz),  however,  the  fluctuations 
become  larger  and  in  some  cases  even  erratic.  This  is  from  system  noise  because  of  the 
low  signal-to-noise  ratio  at  the  edges. 
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At  250  kHz,  the  fast  and  slow  wave  impedanees  are: 

Zj  =(3.072  ±0.097)  10'^  Pa- 

Z2  =(6.2±8.8)-10'^ 

As  ean  be  seen,  the  impedanee  of  the  slow  wave  is  eonsistent  with  zero.  This  plot  does 
not  support  the  existence  of  the  Biot  waves. 
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4.2  WAVE  SPEEDS  AND  MODAE  MASSES 

The  fast  wave  number  is  ealeulated  from  eqs.  (3-2)  and  (3-6): 
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Beeause  the  natural  logarithm  of  a  eomplex  number  is  the  logarithm  of  the  absolute  value 
plus  i  times  the  angle  of  the  complex  number,  and  the  wave  phase  comes  from  the  imagi¬ 
nary  part  of  i  times  the  wave  number,  the  real  part  of  the  fast  wave  number  is: 


=  Re(^i  )=  ^arg(R3^,.,_2«r ), 


(4-8) 


where  the  thickness  is  4  =  1.008  in  or  0.02560  m.  Elsing  MatLAB,  k\  (as  a  function  of 
frequency)  is  “unwrapped”  and  adjusted  so  that  its  intercept  is  zero.  The  fast  wave  speed 
is: 


Cl  = 


(4-9) 


which  is  shown  in  Eig.  4.3. 

The  modal  mass  (effective  density)  of  the  fast  wave  is: 

A=  — •  (4-10) 

Cl 

The  density  of  the  water  filled  sand  is: 

P  =  ^  +  P^  (4-11) 

and  is  plotted  in  Eig.  4.2.  In  eq.  (4-1 1),  is  the  buoyant  weight  of  the  sand,  V  is  the 
total  volume  of  the  water  filled  sand,  and  pw  is  the  density  of  water.  The  diameter  and 
thickness  values  are  averages  at  100  different  locations  measured  to  the  nearest  0.001  in. 
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=14.462  +  0.004  kg 


V  =  71 - ,  D  =  diameter,  =  thickness  of  sand 

4 

Z)  =  3 1 .992  ±  0.001  in  =  0.8 1260  ±  0.00003  m 
=  1 .008  ±  0.001  in  =  0.02560  ±  0.00003  m 

V  =  810.3  ±  0.3  in^  =  0.013278  ±  0.000005 


998.016 


at  20.9°  C. 


The  density  and  speed  of  the  slow  wave  are; 

Pi  =  P~  P\ 


(4-12) 


Cl  —  •  (4-13) 

Pi 

The  magnitudes  of  the  modal  densities  computed  from  eqs.  (4-10)  and  (4-12)  and 
their  standard  deviations  are  shown  in  Fig.  4.2.  The  wave  speeds,  eqs.  (4-9)  and  (4-13) 
and  standard  deviations  are  shown  in  Fig.  4.3.  Porosity  and  sand  density,  which  will  be 
used  later,  are: 

(w) 

fYl  —  fTL^  ' 

R  —  \ 

•  X  7-  _ 


Mg  =mass  of  sand  =  23.076 ±0.004  kg 
=  0.3500  ±0.0002 


(4-17) 


'W"''S  _ 


(w) 


=  2673±2^. 


The  sand  density  compares  to  2650  kg/m^  commonly  used  for  silica. 


(4-18) 


48 


The  fast  and  slow  effective  density  numbers  in  eq.  (4-14)  are  at  250  kHz.  As  can  been 
seen,  the  density  of  the  fast  wave  is  89%  +  3%  of  the  total  density.  The  difference  be¬ 
tween  the  total  density  and  the  fast  effective  density  is  3.9  standard  deviations.  This  fact 
indicates  the  existence  of  Biot  waves  in  water-saturated  sand.  This  compares  favorably 
with  Williams  [42]  who  predicts  a  fast  effective  density  that  is  85%  of  total  density  (see 
Table  5.2). 


1.5  2  2.5  3  3.5 
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Figure  4.3 

Fast  and  Slow  Wave 
Speeds 


At  250  kHz  the  fast  and  slow  wave  speeds  are; 

q  =1654±6m/ 


170  ±  560  m/. 


(4-15) 


This  is  the  first  time  the  speed  of  the  slow  wave  has  been  bounded  in  sand.  The  speed  of 
the  slow  wave  is  consistent  with  zero.  But  this  is  to  be  expected  since  Z2  is  consistent 
with  zero. 

The  fast  wave  to  water  speed  ratio  is: 

q  _1654 


1485 
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The  plus  is  used  in  front  of  the  radical  because  the  answer  for  a  negative  radical 
leads  to  an  inconsistency  in  the  mode  shape  calculation.  In  the  next  step,  mode  shapes 
are  calculated  and  a  negative  radical  produces  a  fast  wave  mode  shape  of  one,  i.e.,  B\IA\ 
=  1.0  +  0.5.  Substitute  BilAi  =  1  back  into  eq.  (2-46)  the  modal  mass  becomes  2087 
rather  than  1850.  Effective  pressures  and  their  standard  deviations  are  shown  in  Fig.  4.4. 
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(4-22) 


The  effective  pressures  at  250  kHz  are; 

Pi  =1.044  ±0.046 
Pj  =-0.044  ±0.046 

The  effective  pressure  of  the  slow  wave  is  consistent  with  zero.  The  measurement  of  ef¬ 
fective  pressures  neither  supports  nor  denies  the  existence  of  Biot  waves. 


4.4  A  SECOND  LOOK  AT  THE  PROPERTIES  OF  THE  FAST  WAVE 

The  effective  density  of  the  fast  wave  being  smaller  than  the  total  density  of  the 
water-saturated  sand  is  the  only  indication  of  the  existence  of  Biot  waves.  The  two  other 
measured  parameters  that  might  confirm  Biot  waves,  the  slow  wave  impedance  and  slow 
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wave  effective  pressure,  are  zero  within  the  standard  deviation.  The  fast  density  is  calcu¬ 
lated  from  impedance  and  speed.  The  importance  of  these  parameters  (fast  wave  imped¬ 
ance  and  fast  wave  speed)  justifies  another  look. 

Fig.  3.4  shows  the  reflection  from  the  second  surface,  which  is  acrylic-air  (r2a,>), 

at  0.92  10”^  sec.  But  the  reflection  from  the  second  surface  in  Fig.  3.2,  acrylic-sand 
{r2s),  at  the  same  time  is  practically  non-existent.  This  implies  the  total  impedance  of  the 
sand  is  equal  to  the  impedance  of  the  acrylic.  In  a  closed  pore  situation,  the  total  imped¬ 
ance  of  the  sand  is  equal  to  the  impedance  of  the  fast  wave  plus  the  impedance  of  the 
slow  wave.  Therefore  the  impedance  of  the  fast  wave  should  be  slightly  less  than  the 

impedance  of  the  acrylic.  The  measured  impedance  of  the  acrylic  is  3.17  -10®  Pa  ■  s/m, 
see  eq.  (A-6).  Ref  [27]  measures  the  speed  of  sound  in  cast  acrylic  at  2705  m/s  and  the 
density  to  be  1181  kg/m^,  which  makes  the  impedance  3.19  10^.  The  measured  imped¬ 
ance  of  the  fast  wave  in  this  report  is  smaller  at  3.07  10^  Pa-s/m,  which  is  consistent 
with  expectations. 

The  speed  of  the  fast  wave  in  Fig.  4.3  was  found  to  be  c\  =  1654  m/s  by  phase 
measurements  in  the  frequency  domain.  In  the  time  domain  the  reflection  from  the  third 
surface  in  Fig.  3.6,  which  is  a  sand-Ethafoam  interface  {tse),  arrives  at  a  time  of 

1.23-10”^  sec.  This  makes  the  time  difference  {r^E  minus  r2air,  f'2air  from  above)  of 

0.31- 10“^  sec.  The  distance  traveled  is  2  times  0.02560  m,  which  gives  a  wave  speed  of 
1650  m/s,  consistent  with  the  frequency  domain  result. 

It  can  be  concluded  that  the  measured  values  for  the  fast  wave  impedance  and 
speed  of  the  fast  wave  are  reasonable  and  consistent  with  expectations. 
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4.5  MODE  SHAPES 

Eig.  4.5  is  a  model  of  a  closed  pore-sediment  interface  for  the  first  mode.  If  the 
incident  and  reflected  waves  on  the  interface  have  a  combined  amplitude  of  one  (which 
makes  the  first  mode  displacement  u\  =  l-exp(/(tt)t  -  k\z))  ),  then  from  eqs.  (2-26),  (2-27), 
(2-48),  and  (2-49): 

A,{\-P)+B,P  =  P,  ,  (4-23) 

and  likewise: 

A^{\-P)+B^fi  =  P^  .  (4-24) 

The  participation  factors.  Pi  and  P2 ,  are  independent  of  boundary  conditions  but  they 
just  happen  to  be  (luckily)  numerically  equal  to  the  open  pore  effective  forces.  Pi  and  P2, 
eq.  (4-22).  Getting  back  to  the  closed  pore  interface,  Ai  +  A2  =  Bi  +  B2,  which  means 
there  is  no  relative  movement  between  sand  and  fluid  in  the  sediment  making  viscous 
forces  and  tortuosity  of  the  fast  wave  cancel  viscous  forces  and  tortuosity  of  the  slow 
wave  at  the  surface.  Therefore  the  mass  that  moves  is: 

(4-25) 

A2{l-fi)p,+B2fiPf=P2-  (4-26) 
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tor  East  W  ave 

Solving  for  the  mode  shapes,  i.e.,  the  sediment  fluid  displacement  divided  by 
sand  displacement: 
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The  fast  wave  mode  shape  is  3.2  standard  deviations  from  one  and  helps  support  the  ex- 
istenee  of  Biot  waves. 

The  fast  wave  mode  shape  in  Fig.  4.6  is  an  important  measurement  in  this  disser¬ 
tation.  This  result  means  that  the  water  in  the  sediment  oscillates  with  an  amplitude  (and 
velocity  and  acceleration)  2.2  times  the  amplitude  of  the  sand  during  the  passing  of  a 
sound  wave.  A  difference  in  vibration  amplitudes  between  sand  and  fluid  is  one  of  the 
key  properties  of  a  Biot  wave.  This  fact  gives  a  clearer  understanding  of  the  physics  of 
sound  waves  propagating  through  sandy  sediment  more  than  any  other  parameter  meas¬ 
ured  in  this  dissertation. 


4.6  TORTUOSITY  AND  VISCOSITY 


The  tortuosity  and  viscosity  terms  always  appear  together  and  cannot  be  solved 
for  separately.  Eq.  (2-14)  can  be  rewritten: 


~P  Q 

Q  R. 

Ppf-a 

Ppf  -a  a 


\b\ 


(4-31) 


Where  the  “mass  factor”  is  (units  of  mass  or  density): 

a  =  aPp f  -i  — 
CO 


and  will  be  solved  for  in  the  next  section. 


(4-32) 


4.7  MODULI 

From  eq.  (4-31)  P,  Q,  R,  and  a  can  be  solved  for: 
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Figure  4.8 

Mass  Factor 

At  250  kHz  the  mass  factor  is: 

a  =  340  ±110  kg.  (4-37) 

Using  the  known  values  of  porosity  and  water  density  eq.  (4-32)  becomes: 

a  =  M9a-i—.  (4-37a) 

CO 

If  ct=  1  and  F  =  0,  i.e.,  no  tortuosity  or  viscous  effects,  then  a  in  eq.  (4-37a)  is  0.08 
standard  deviations  from  eq.  (4-37).  This  agrees  with  an  experiment  previously  refer¬ 
enced,  [33],  which  concluded  “attenuation  shows  no  correlation  with  the  viscosity”  of  the 
fluid  in  unconsolidated  sand. 
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Note  that  up  to  now  all  graphed  parameters  are  measured.  Starting  with  the  pres¬ 
sure  refleetion  pulses  and  ending  with  the  generalized  moduli  and  mass  faetor,  all  math 
and  theory  (assuming  the  equations  of  motion  and  boundary  eonditions  are  true)  are  rig¬ 
orous.  But  the  following  sand,  frame,  and  shear  moduli  formulas  are  “semiphenome- 
nologieal”  [20]  and  not  rigorous. 

From  eqs.  (2-5)  and  (2-4)  the  sand  and  sand  frame  moduli  are; 


QK, 


j3K^-R 


Kb  = 


Ql 

R 


Ks- 


(4-38) 


(4-39) 


The  modulus  of  water,  ifw,  is  =  (l485.l)^  (998.01)  =  2.2012  ■  10^  Pa  at  20.9°  C. 


From  eq.  (2-3)  the  shear  modulus  is: 


(4-40) 


Moduli  at  250  kHz  are; 


=(7.3 ±3.1)- 10^  Pa 

(4-41) 

/V  =  (9.0  ±7.5)  10^  Pa 

(4-42) 

Kf,  =(-5.6 ±7.3)- 10^  Pa 

(4-43) 

These  ealeulations  do  not  appear  eorreet.  is  9  standard  deviations  lower  than 
the  reeognized  value  for  sand  modulus  (about  36  GPa  [32]).  It  is  interesting  that  Choti- 
ros  [13]  also  ealeulated  a  sand  modulus  of  Ks  =  7  GPa  using  a  numerieal  inversion  teeh- 
nique  from  measured  parameters  (he  eoneluded  7  GPa  was  low).  The  formulas  for  the 
generalized  moduli,  eqs.  (2-3),  (2-4),  and  (2-5),  assume  the  bonding  of  the  grains  to  be  of 
the  same  material  as  the  grains,  ref  [20]  p.  558.  This  assumption  may  not  be  eorreet  for 
uneonsolidated  sand.  The  generalized  moduli  formulas  appear  eorreet  when  the  grains 
are  bonded  [22],  but  the  “bonding”  of  sand  grains  may  inelude  sliding  frietion  and/or  a 
fluid-reinforoed  frame  [7,  11].  Finding  good  formulas  for  intermediate  moduli  of  sandy 
sediments  is  a  task  for  future  work. 


4.8  OPEN  PORE  PRESSURE  REEEECTION  COEEEICIENTS 
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The  conventional  reflection  coefficient  using  the  density  of  the  water  saturated 
sand  and  the  speed  of  sound  in  sand; 

R  =  =  0.40  .  (4-47) 

Pci+Z^ 

These  reflection  coefficients  compare  to  a  reflection  experiment  Chotiros  et  al. 
[12]  did  on  sandy  sediment.  His  reflection  loss  was  1 1  +  2  dB  (R^v  =  0.28  +  0.08),  to  be 
compared  to  a  computed  conventional  reflection  loss  of  8  +  1  dB  (7?  =  0.40  +  0.05). 


61 


Chapter  5: 

DISCUSSION 

This  was  an  experimental  dissertation  to  help  understand  the  physics  of  sound 
propagating  through  sandy  sediments.  The  measurements,  which  were  obtained  through 
a  reflection  ratio  technique,  include  acoustic  impedances,  effective  densities,  waves 
speeds,  effective  pressures,  mode  shapes,  intermediate  moduli,  and  pressure  reflection 
coefficients.  A  lumped  parameter  model  was  developed  to  help  aid  in  the  interpretation 
of  the  measurements.  The  results  show  that  the  effective  density  of  the  fast  wave  is  less 
than  the  total  density  of  the  water-saturated  sand.  This  fact  points  to  the  existence  of 
Biot  waves  in  unconsolidated  sand. 

Three  characteristics  of  a  fast  Biot  wave  are  that  it  produces  different  water/sand 
oscillation  amplitudes,  that  it  does  not  use  the  full  inertial  energy  in  the  sediment,  and  it 
has  a  different  phase  speed  from  its  slower  Biot  wave  even  though  both  are  compres- 
sional  waves  traveling  through  the  same  medium. 

5.1  COMPARISONS 

In  the  section  that  follows  selected  parameters  are  calculated  from  Biot  theory  at 
250  kHz  using  published  input  parameters  and  compared  to  parameters  measured  from 
this  experiment.  Williams’  parameters,  ref  [42]  in  Table  5.1  are  used  to  calculate  se¬ 
lected  parameters  that  are  compared  to  parameters  measured  in  this  experiment.  The 
comparisons  are  tabulated  in  Table  5.2.  The  deviations  in  Table  5.2  are  the  number  of 
standard  deviations  the  measured  differs  from  theoretical.  No  deviations  were  listed  for 
mechanically  measured  parameters.  As  previously  mentioned  the  sand  modulus  is  over  9 
standard  deviations  lower  than  expected,  presumably  due  to  incorrect  intermediate 
moduli  formulas.  The  result  of  this  is  the  large  deviation  of  the  fast  wave  speed. 

The  fast  effective  mass  using  Williams’  parameters  and  lumped  parameter  model 
is  (Table  5.2): 

Pi  =1692.2-/6.5  .  (5-1) 

/  m 
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This  compares  with  the  effeetive  mass  using  Williams’  parameters  and  effeetive  density 
method  [42]; 

=1674.2-/1.2  V3  (5-2) 

/  m 

The  sign  on  the  imaginary  part  is  negative  beeause  of  the  forms  of  the  solutions  assumed, 
eqs.  (2-21)  and  (2-22).  The  differenee  between  eqs.  (5-1)  and  (5-2)  is  due  to  the  ap¬ 
proximations  used  in  [42]. 


Table  5.1 

Williams’  Parameters  rei 

f.  [42] 

Physieal  Constants 

Symbols 

Units 

Values 

Permeability,  see  note  1 . 

k 

m^ 

le-10 

Porosity 

P 

— 

0.40 

Density  of  fluid 

Pf 

kg/m"* 

1000 

Density  of  sand  grains 

Pr 

kg/m"* 

2650 

Bulk  modulus  of  grains 

Ks 

Pa 

3.6el0 

Bulk  modulus  of  fluid 

K„ 

Pa 

2.25e9 

Bulk  modulus  of  frame,  see  note  2. 

Kh 

Pa 

4.4e7  +  i2.0e6 

Fluid  viseosity 

V 

Ns/m"^ 

le-3 

Shear  modulus,  see  note  2. 

N 

Pa 

2.61e7  +  il.24e6 

Virtual  mass  for  liquid  (tortuosity) 

a 

— 

1.25 

Pore  size  (radius) 

a 

m 

jSak 

i  P 

Notes: 

1 .  Wave  number  and  permeability  both  use  k. 

2.  The  signs  on  the  imaginary  parts  are  positive  beeause  of  the  forms  of  the  solu¬ 
tions  assumed,  eqs.  (2-21)  and  (2-22). 
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Table  5.2 

Comparisons  at  250  kHz 

Parameter 

Measured 

Williams’ 

Deviation 

(Symbols) 

Parameters  [42] 

See  note  1 

Fast  Impedance 

(Zi) 

(3.072  +  0.097)e6 

(2.967  -  z0.005)e6 

-1.1 

Slow  Impedance 

(Z2) 

(6.2  ±  8.8)e4 

(4.66  +  z0.26)e4 

0.2 

Total  Density 

(P) 

2087 

1990 

— 

Fast  Effective  Density 

(A) 

1860  +  58 
89%±3%  of  total 

1692.2  -  i6.5 

85%  of  total 

-1.3 

dev.  of  % 

Slow  Effective  Density 

(A) 

227  +  58 

1 1%±3%  of  total 

297.5  + 16.5 

15%  of  total 

1.3 

dev.  of  % 

East  Wave  Speed 

(ci) 

1654  +  6 

1753.4 +  z3.8 

see  note  2 

Slow  Wave  Speed 

(C2) 

170  +  560 

156.7  +  z5.4 

-0.02 

East  Effective  Eorce 

{Pi) 

1.044  ±0.046 

1.0029 +  Z0.0002 

-0.9 

Slow  Effective  Eorce 

{Pi) 

-0.044  +  0.046 

-0.0029  -  zO.0002 

0.9 

East  Mode  Shape  {Bi/Ai) 

2.20  +  0.37 

2.100 +  Z0.036 

-0.3 

Slow  Mode  Shape  (B2/A2) 

-2.37  ±0.47 

-1.5235 -zO.0008 

1.8 

Intermediate  Modulus 

{P) 

(2.69  ±  0.45)e9 

(1.923  +  z0.003)e9 

-1.7 

Intermediate  Modulus 

(0 

(9.7±  l.l)e8 

(1.2320 -z0.0001)e9 

2.4 

Intermediate  Modulus 

{R) 

(4.8±  1.0)e8 

8.2300e8  +  z7e3 

3.4 

Mass  Eactor 

{a) 

340  ±  no 

511-/12 

1.6 

Sand  Modulus 

(^s) 

(7.3±3.1)e9 

3.6el0 

9.3 

Erame  Modulus 

(^b) 

(-5.6±7.3)e8 

(4.4  +  /0.2)e7 

0.8 

Shear  Modulus 

(AO 

(9.0±7.5)e8 

(2.61  +  /0.12)e7 

-1.2 

Porosity 

{P) 

0.35 

0.40 

— 

Sand  Density 

(A) 

2673 

2650 

— 

Water  Speed 

(S^w) 

1485 

1500 

— 

Note  1:  deviation  =  (Williams’  minus  measured)/standard  deviation 

Note  2:  A  better  comparison  might  be  the  fast  speed  to  water  speed  ratios. 

Por  Meas- 

ured  this  is  1.11.  Prom  Williams’  parameters  it  is  1.17.  A  difference  of  5%. 
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Table  5.3  shows  the  units  of  the  parameters. 


Table  5.3 

Parameter  Units 

Parameter 

Symbol 

Continuous  Units 

Eumped  Units 

Fast  Impedance 

Zi 

Pa-s/m 

Ns/m 

Slow  Impedance 

Z2 

Pa-s/m 

Ns/m 

Total  Mass 

P 

kg/m"* 

kg 

Fast  Effective  Mass 

pi 

kg/m"* 

kg 

Slow  Effective  Mass 

pi 

kg/m"* 

kg 

East  Wave  Speed 

Cl 

m/s 

m/s 

Slow  Wave  Speed 

C2 

m/s 

m/s 

East  Effective  Eorce 

Pi 

Pa/Pa 

N/N 

Slow  Effective  Eorce 

Pi 

Pa/Pa 

N/N 

East  Mode  Shape 

Bi/Ai 

m/m 

m/m 

Slow  Mode  Shape 

B2/A2 

m/m 

m/m 

Intermediate  Modulus 

P 

Pa 

N/m 

Intermediate  Modulus 

Q 

Pa 

N/m 

Intermediate  Modulus 

R 

Pa 

N/m 

Mass  Eactor 

a 

kg/m"* 

kg 

Sand  Modulus 

Ks 

Pa 

N/m 

Erame  Modulus 

Kb 

Pa 

N/m 

Shear  Modulus 

N 

Pa 

N/m 

Porosity 

P 

m'/m' 

m/m 

Sand  Density 

A 

kg/m"* 

kg 
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5.2  SIMPLIFIED  EQUATIONS  OE  MOTION 

At  higher  frequencies  (250  kHz)  for  water-saturated  sand  at  normal  incidence  the 
equations  of  motion,  eqs.  (2-1)  and  (2-2)  reduce  to  (no  viscous  or  tortuosity  effects): 

d^u 


P^-^  +  Q^-{l-P)PsU,=0 

02 


Q 


dh 


+  R- 


V 


PpfUf  =  0. 


(5-3) 


(5-4) 


dz^  dz^ 

These  simplified  equations  only  use  seven  of  the  13  Biot  parameters  listed  in  Ta¬ 
ble  5. 1 :  (5,  pi,  fh,  As,  Aw,  and  Ab  (complex  counts  as  two). 


5.3  EUTUREWORK 

Euture  work  should  include  (1)  adding  shear  displacements  to  the  lumped  pa¬ 
rameter  model  (allowing  for  angle  of  incidence),  (2)  deriving  good  intermediate  moduli 
formulas,  (3)  using  the  reflection  ratio  technique  on  a  different  unconsolidated  medium 
like  glass  beads.  This  may  help  to  derive  good  moduli  formulas.  (4)  The  development 
of  a  “calibrated”  piece  of  material,  e.g.,  cast  acrylic,  that  can  be  taken  into  the  field  and 
put  on  top  of  the  sand.  With  a  transducer  and  receiver  mounted  a  fixed  distance  from  the 
acrylic,  coherently  averaged  fast  wave  impedances  can  be  measured  (may  have  to  as¬ 
sume  slow  wave  impedances  are  small).  That  way,  together  with  the  speed  of  sound,  ef¬ 
fective  densities  of  the  sediment  can  be  measured  in  the  field. 


5.4  CONCEUSIONS 

This  dissertation  is  a  data  point  in  favor  of  Biot  waves  existing  in  sandy  sedi¬ 
ments  because  (1)  the  fast  wave  effective  density  is  less  than  the  total  density,  (2)  the  fast 
wave  mode  shape  is  not  one,  and  (3)  the  coupling  modulus,  Q,  is  not  zero.  However, 
other  measured  parameters,  e.g.,  slow  wave  impedance  and  speed,  slow  wave  effective 
pressure,  tortuosity  and  viscous  effects,  are  consistent  with  zero  and  do  not  support  the 
existence  of  Biot  waves.  Even  the  Biot  open  pore  pressure  reflection  coefficient  is  con¬ 
sistent  with  the  conventional  pressure  reflection  coefficient. 
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The  largest  eontribution  toward  understanding  the  physies  of  sound  waves  in 
sandy  sediments  is  the  measurement  of  mode  shapes,  whieh  give  relative  oseillation  am¬ 
plitudes  between  the  saturating  fluid  and  solid  partieles.  The  knowledge  that  the  sand 
and  fluid  vibrate  with  different  amplitudes  helps  explain  how  a  sound  wave  propagates 
without  using  the  full  mass  of  the  medium.  But  the  largest  praetical  contribution  is  the 
measurement  of  the  effective  density  of  the  fast  wave.  Effective  densities  can  be  used  for 
more  accurate  predictions  of  scattering  [44]. 

This  experiment  supports  the  conclusions  that  subcritical  penetration  in  sandy 
sediment  is  from  scattering  and  not  from  the  Biot  slow  waves  [17,  26,  34].  This  is  true 
for  two  reasons.  The  speed  of  the  slow  wave  is  essentially  zero  and  the  incident  pressure 
wave  does  not  excite  the  slow  wave  (participation  factor  is  essentially  zero).  In  other 
words,  no  energy  goes  into  a  slow  wave  and  if  it  did  it  would  not  propagate. 
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Appendix  A: 

IMPEDANCES  OE  ACRYEIC  AND  ETHAEOAM 

A.l  IMPEDANCE  OE  ACRYEIC 

The  measurements  and  ealculations  for  the  impedanee  of  the  acrylic  were  accom¬ 
plished  using  ratios  of  transfer  functions.  The  acrylic  was  purchased  from  Port  Plastics, 
Tukwila  WA  under  the  product  name  of  “Clear  Acrylic.”  The  transfer  functions  were  the 
ratios  of  reflections  shown  in  Pigs.  A.l  and  A.2.  In  these  figures  the  transducer  and  re¬ 
ceiver  are  in  the  water  below  the  acrylic. 


The  received  pressure  from  the  transmitter,  reflected  off  the  water-acrylic  inter¬ 
face,  and  back  to  the  receiver  is: 

Ria=^^ — (A-1) 
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A  is  the  complex  pressure  amplitude  a  unit  distance  from  the  acoustic  center,  di  is  the 
distance  from  the  acoustic  center,  Za  is  the  impedance  of  the  acrylic,  is  the  impedance 
of  water,  k„  is  the  wave  number  of  water,  and  d„  is  the  distance  the  sound  travels  through 
water  to  the  acrylic  surface  and  back.  Due  to  the  mechanics  of  the  transducers  d\  may 
not  equal  d„.  Riair  is: 

Riair  - ^exp(-zk^flf^)exp(-z2k^tj.  (A-2) 

Where  d2  is  the  spreading  loss  factor,  ka  is  the  wave  number  of  the  acrylic,  and  ta  is  the 
thickness  of  the  acrylic.  The  reflection  coefficient  i?2w  is: 


^iw  - 


A  2Z. 


7  -  7 


27., 


exp(-  zk^  )exp(-  ilkj^ ) . 


The  ratio  of  transfer  functions  for  the  acrylic-air  surface  is: 


(^2^)2 
i^la  ^2 
i^2air  )i 
{Rla  )l 


^ 7  -7  ^ 

\^w  J 


(A-3) 


(A-4) 


The  subscripts  1  and  2  represent  acoustic  pulses,  in  which  the  transducer-generated  am¬ 
plitude,  A,  and  the  distance  through  the  water,  dw,  may  be  slightly  different  from  one 
setup  to  the  next.  Also  7?2w  and  Rimr  are  divided  by  7?ia  to  make  the  reflection  from  the 
acrylic  surface  the  t  =  0  trigger.  Solving  eq.  (A-4)  for  7^: 


(i  +  r)z.. 


(A-5) 


The  impedance  and  standard  deviation  at  100  locations  of  the  4”  acrylic,  Za,  are  shown  in 
Fig.  A.3. 


1  1.5  2  2.5  3  3.5  4 


frequency  (Hz) 

Figure  A.  3 
Impedance  of 
Acrylic 

For  the  acrylic,  at  250  kHz,  the  impedance  is: 

=(3.17±0.05)-10'^  Pa-s/m.  (A-6) 

As  a  comparison,  ref.  [27]  measures  the  speed  of  sound  in  cast  acrylic  at  2705  m/s  and 

3  6 

the  density  to  be  1181  kg/m  ,  which  makes  the  impedance  3.1910  . 
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A.2  IMPEDANCE  OE  ETHAEOAM 

Ethafoam  is  similar  to  Styrofoam.  Ethafoam  is  a  little  denser  but  is  easier  to 
work  with  and  glue.  The  transfer  funetions  (ratios  of  refleotions)  are  shown  in  Eigs.  A.l 
and  A.4. 


The  ratio  of  transfer  funetions  for  the  aerylie-Ethafoam  surfaee  is: 


^lE 

^\a 

(^Ze  -Za 

^lair 

^\a 

(A-7) 


Where  Ze  is  the  impedanee  of  the  Ethafoam  and  is  the  impedanee  of  the  4”  aerylie. 
Solving  eq.  (A-7)  for  Zg. 

Ze  =  ~  (A-8) 


The  impedanee  of  the  Ethafoam,  Ze,  and  its  standard  deviation  are  shown  with  in  Eig. 
A.5. 


Z£  is  0.5  standard  deviations  from  zero.  A  Ze  of  zero  is  used  in  the  body  of  this  report. 
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Appendix  B: 
NUMERICAL  CHECKS 


The  lumped  parameter  method  has  two  areas  of  new  theory,  where  the  final  prod¬ 
uets  are  pressure  eoeffieients  and  mode  shapes;  these  two  parameters  are  ealeulated  from 
eontinuous  and  lumped  parameter  theory  using  Williams’  parameters.  Even  though  this 
is  not  a  suffieient  proof  that  the  lumped  and  eontinuous  parameter  theories  are  equivalent, 
it  shows  the  two  methods  are  equivalent  in  the  range  of  values  measured  in  this  experi¬ 
ment. 

Pressure  eoeffieients  eomputed  using  Johnson  and  Fiona’s  method  [20]  are  ealeu¬ 
lated  below  and  eompared  to  pressure  eoeffieients  from  the  lumped  parameter  method. 
Table  2.1.  Also  modes  shapes  from  Biot  theory,  eqs.  (2-18)  and  (2-19),  are  ealeulated 
and  eompared  to  modes  shapes  from  the  lumped  parameter  method,  eqs.  (4-27)  and  (4- 
28).  All  pressure  eoeffieients  and  mode  shapes  are  ealeulated  using  the  input  parameters 
of  Table  5.1. 

As  the  ineident  pressure  from  the  water  aets  upon  the  sediment,  the  foree  in  the 
sand  per  unit  area  of  sediment  equals  the  foree  on  the  sand  from  the  ineident  and  re- 
fieeted  wave  per  unit  area  of  sediment: 


+  0 


dz 


i^-mK 


^Uj 

dz 


+  K^ 


(B-1) 


And  likewise,  the  foree  in  the  fluid  (the  fluid  in  the  sediment)  per  unit  area  of  sediment 
equals  the  foree  on  the  fluid  from  the  ineident  and  refieeted  wave  per  unit  area  of  sedi¬ 
ment: 


'Z  uZ  V  uZ  uZ 


(B-2) 


Another  boundary  eondition  is  partiele  displaeement  below  equals  partiele  displaeement 
above: 

—  I3)u^  +  puj-  =Ui  +u^.  (B-3) 

The  forms  of  the  solutions  are: 
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=Ai  exp(i((i)t - kiz})+  A2  exp(i((i)t - k2z}} 
u  j-  =  Bi  exp(/(ftjf  —  A:jz))+ ^2  Qxp{i{cot  —  ki^)) 
Uj  =  Ai  Qxp{i{cot - k^z)) 

=  A^  exp(/(^yf  +  A:^z)). 


The  five  unknowns,  Ai,  Aj,  B\,  B2,  and  A^,  are  solved  in  matrix  form,  assuming  At  =  1. 
The  five  equations  are  (B-1),  (B-2),  (B-3),  (2-18),  and  (2-19),  yielding 


'A 

~  Pk^ 

Pk2 

Qkx 

Qk2 

(1-AaA" 

-1 

'(1-AA^w' 

A 

Qki 

Qk2 

Rk^ 

Rk2 

PK^k^ 

Pkyv^yv 

<B, 

.  = 

P 

P 

-1 

* 

1 

Bi 

M41 

0 

M43 

0 

0 

0 

A, 

0 

M52 

0 

M54 

0 

0 

where 


M41  =  ki  P  -  aP' Pii  -icoF 
Mat,  =k]Q-co^P]2  +icoF 
M52  =k2P  -  0)^ pii  -icoF 
M54  =  A:|g  -  cd^Pi2  +  i(oF. 

The  equation  (B-3)  was  multiplied  by  10e5  to  keep  the  inversion  in  eq.  (B-5) 
from  becoming  singular. 

The  pressure  reflection  coefflcient  is 


The  transmitted  pressures  coefflcients  of  the  fast  and  slow  waves  are: 


Pi 


h 

k  K 


[{P  +  Q)A,+{Q  +  R)B,] 


Pi 


-^[{p  +  Q)A2+{Q  +  R)B2l 

W^W 


(B-7) 


(B-8) 

(B-9) 


Pressure  coefflcient  comparisons  are  in  Table  B.l. 
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Table  B.l 

Comparisons  of  Caleulated 

Pressure  Coeffieients 

PrIPi 

eq.  (B-7) 

0.32564926491560  -  iO.00091523022991 

Rww 

eq.  (2-57) 

0.32564926491560  -  iO.00091523022991 

GxIPi 

eq.  (B-8) 

1.32943502471988  -  iO.00069866204680 

Twi 

eq.  (2-59) 

1.32943502471988  -  iO.00069866204680 

OltPi 

eq.  (B-9) 

-0.00378575980428  -  iO.00021656818311 

Tw2 

eq.  (2-60) 

-0.00378575980428  -  iO.00021656818311 

Mode  shape  eomparisons  (eontinuous  eqs.  (2-18),  (2-19)  and  lumped  eqs.  (4-27), 
(4-28))  using  input  parameters  in  Table  5.1  are  in  Table  B.2. 

Onee  again,  even  though  these  agreements  in  Tables  B.l  and  B.2  are  neeessary, 
they  are  not  a  suffieient  proof  that  the  lumped  parameter  method  is  equal  to  the  eontinu¬ 
ous  parameter  method.  It  does  show  the  two  methods  are  numerieally  equivalent  in  the 
parameter  range  of  interest. 


Table  B.2 

Comparisons  of  Caleulated 

Mode  Shapes 

By! Ax 

eq.  (2-18) 

2.10033988149722  +  iO.03549735086698 

BxlAx 

eq.  (4-27) 

2.10033988149741  +  iO.03549735086701 

B2IA2 

eq.  (2-19) 

-1.52352553508303  -  iO.00083813014479 

B2IA2 

eq.  (4-28) 

-1.52352553508303  -  iO.00083813014479 
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Appendix  C: 

ERROR  ANALYSIS 

All  acoustically  measured  parameters  were  measured  in  100  different  locations. 
The  bands  or  tolerances  used  in  the  plots  in  the  body  of  this  report  are  standard  devia¬ 
tions.  The  following  calculates  the  errors  at  the  center-band  frequency  (250  kHz)  and 
compares  them  to  the  standard  deviations. 

Three  types  of  error  are  calculated,  random,  bias,  and  “absolute,”  on  Z\  and  Z2. 
Random  error  of  a  random  quantity  is  the  standard  deviation  divided  by  the  square  root 
of  the  number  of  measurements.  So  from  eq.  (4-5),  the  random  error  on  the  fast  and  slow 
wave  impedances  are  9700  and  8800  Pa- s/m  respectively. 

The  residual  ringing  of  the  transmitting  transducer  caused  a  bias  error,  see  Ap¬ 
pendix  G.  Fig.  C.l  is  the  pulse  reflected  from  the  water-air  interface  along  with  the  en¬ 
velope  of  the  pulse  for  acrylic-air,  see  Fig.  3.4.  The  envelope  is  placed  at  the  same 
time  the  second  pulse  from  the  top  acrylic  surface  would  return.  The  value  under  the  en¬ 
velope  of  interest  is  used  for  the  bias  errors  for  7?3e  3air  and  7?2s_2air  in  eqs.  (4-3)  and  (4-4). 
The  bias  errors  of  Z\  and  Z2  are  calculated  using  the  partial  derivatives  of  eqs.  (4-3)  and 
(4-4)  with  respect  to  R3E  3air  and  7?2s_2air  and  the  bias  errors  of  7?3e jair  and  R2s_2air-  The 
bias  errors  of  Z\  and  Z2  are  17,000  and  7500  Pa- s/m  respectively. 

The  “absolute”  error  arises  from  the  absolute  values  used  instead  of  the  complex 
on  the  reflections  in  eqs.  (4-1)  and  (4-2).  To  obtain  absolute  errors,  William’s  parame¬ 
ters  are  used  to  calculate  complex  Z\  and  Z2,  as  in  Table  5.2.  The  complex  reflections  are 
calculated  from  eqs.  (4-1)  and  (4-2),  made  into  absolutes,  and  absolute  Z\  and  Z2‘s  are 
calculated  with  eqs.  (4-3)  and  (4-4).  The  “absolute”  errors  are  found  to  be  16,000  and 
21,500  Pa- s/m  respectively. 

Assuming  these  errors  are  not  correlated  and  can  be  added  in  quadrature,  the  total 
errors  of  Z\  and  Z2  are  25,000  and  24,000  Pa-s/m  respectively  (relative  error  0.008  and 
0.39),  compared  to  standard  deviations  of  97,000  and  88,000  Pa-s/m  respectively  (rela¬ 
tive  deviation  0.03  and  1.4). 


amplitude  squared  amplitude 
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The  total  errors  are  about  one-third  the  standard  deviations.  Therefore  the  stan¬ 
dard  deviations  are  mostly  a  measure  of  spatial  variability  of  the  sand  sample. 


0  0.5  1  1.5  2 


Figure  C.l 

Pulse  Reflection  from  Water-Air  Interface 
with  Envelope  Used  for  Bias  Error 
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Appendix  D: 

SINGLE  DEGREE-OF-FREEDOM  EUMPED  PARAMETER  MODEE 

The  first  part  of  this  section  shows  that  the  motion  from  a  lumped  mass,  linear 
spring  model  is  identical  to  the  motion  of  a  particle  in  a  continuous  medium  where  there 
are  no  boundary  effects.  The  second  part  does  the  same  for  the  boundary. 


D.l  HOMOGENEOUS  (NO  BOUNDARY  EFFECTS) 

Start  by  taking  a  single  degree-of-freedom  system,  i.e.,  only  one  wave  traveling 
through  the  medium,  and  find  the  homogeneous  solution,  see  Fig.  D.l  below. 


Where  Ki  is  the  linear  stiffness,  pe  is  the  lumped  mass.  The  conversions  of  continuous 
parameters  to  lumped  parameters  are: 


Assume 


Ki  (N/m)  =  ^ 

Az(mj 

/7^(kg)=  /7(kg/m^)xl(m^)xAz(m). 
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(D-1) 
(D-2) 

u_i=  A  Qx^{i(o{t  +  t)) 

mq  =  ^exp(/^yr)  (D-3) 

Ui=  A  QXY){ico{t  -  r)). 

Where  ris  a  time  delay. 

The  motions  of  the  masses  u.\  and  mq  and  their  differences  are  shown  in  Fig.  D.2 
as  a  function  of  time. 


Movements  of  Two  Masses  {u.\  and  mq) 
and  the  Difference  of  Their  Movements 
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The  slope  at  the  zero  crossing  for  w.i  and  m  is  Ao).  The  amplitude  at  a  time  t/2 
from  the  zero  crossing  isAoJtll  when  ris  small.  The  difference  between  uo  and  u.\  is: 


r 

ico 

t+- 

V 

K  2  jy 

Likewise: 


Uq  -u_i  =-iAa>TQxp 

Ui  -Uq  =-iAcuTQxp 


r 

ico 

t — 

V 

'v  2  jy 

And:  (wj  -Uq)-{uq  -  u_i)  =  -Acd^t  Qxp{ici}t). 

The  equation  of  force  on  the  mass  at  mq  is  (exp(zT<tf)  suppressed): 

[(mi  -uq)-{mq  -u_i)]k^  =-Apf(u^. 
Substitute  in  eqs.  (D-1),  (D-2),  and  (D-5): 

_K 

P  ' 

The  speed  of  the  wave  (c)  is  Az/ r,  therefore: 

fY 

ip' 


(D-4) 

(D-5) 

(D-6) 

(D-7) 

(D-8) 


(D-9) 


This  is  the  speed  of  a  wave  in  a  continuous  medium.  The  speed  of  the  wave  is  the 
same  for  a  continuous  medium  as  it  is  for  a  lumped  medium  given  Az  is  much  smaller 
than  the  wavelength. 


D.2  BOUNDARY 

The  pressure  on  the  boundary  (z  =  0)  for  a  continuous  medium  is  the  modulus 
times  the  first  spatial  derivative  (strain). 


du 

p  =  K  — 
dz 


(D-10) 


p  =  K{-iAk  Qxp{iCDt)) 

This  is  assuming  the  movement  of  a  particle  in  the  continuous  medium  is  of  the  form: 

u  =  A  Qxp{i{ciM  -  kz)) . 


(D-11) 
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The  force  on  the  boundary  for  a  lumped  medium  system  is  the  stiffness  times  the 
difference  in  displacement  of  two  points  Az  apart. 


f  =  K({ui  -Mo) 
/  =  -iAkK  exp 


r 

id) 

t — 

V 

'v  2  jy 

(D-12) 


r,  on  0)  ,  1 ^ 

Where;  —  =  —  =  k.  (D-13) 

Az  c 

As  Az  ^  0,  and  therefore  as  0,  the  pressure  and  the  force  on  the  boundary  are 
the  same.  As  pointed  out  earlier  stiffness  (N/m)  and  modulus  (N/m  ),  mass  (kg)  and  den- 
sity  (kg/m  ),  force  (N)  and  pressure  (N/m  ),  viscous  damping  (Ns/m)  and  body  viscous 
damping  (Ns/m"^),  difference  in  displacement  of  two  points  (m)  and  strain  (m/m),  are 
used  interchangeably.  Hopefully  this  does  not  cause  too  much  confusion. 

This  appendix  shows  that  the  continuous  medium  and  lumped  medium  are 
equivalent  for  a  single  degree-of-freedom  system.  This  dissertation  assumes  this  equiva¬ 
lency  holds  for  multiple  degree-of-freedom  systems. 


85 


Appendix  E: 

DYNAMIC  DESIGN  ANAEYSIS  METHOD 


Eigs.  E.I,  E.2,  E.3,  and  E.4  are  eopies  of  selected  pages  of  Belsheim  and 
O’Hara’s  NAVSHIPS  250-423-30,  Shock  Design  of  Shipboard  Equipment,  Dynamic  De¬ 
sign  Analysis  Method  (DDAM),  ref.  [2],  DDAM  analyzes  shock  on  Naval  ships  from 
non-contact  explosions.  Belsheim  and  O’Hara  were  the  first  to  derive  expressions  for 
effective  mass  and  participation  factors.  Modal  mass  is  also  referenced  in  Clough  and 
Penzien  [14]  for  earthquake  analysis  of  structures. 

Rewriting  the  effective  mass  from  DDAM  ([2]  eq.  (7))  and  Clough  and  Penzien, 
p.  559  [14]  to  include  mass  coupling  (for  the  first  mode  only,  second  mode  similar): 


jifUi  ">12 ]|i,  I 

[ij  [^21  «22j[^J 

D  mil  ^12  n 

Yi\  L"^21  ^2i\\Ai 


(E-1) 


Which  is  the  same  as  eq.  (2-46)  if  the  mass  matrix  is: 
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— 

ft) 
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The  DDAM  participation  factor  ([2]  eq.  (6))  is: 


jifKi  ">12]  i, , 

[ij  [m2i  «22jUj 
Bi  ,  "^11  "^12  ,  Bi  , 

Yi\  L'”21  "^22j[)4^ 


(E-2) 


(E-3) 


This  is  different  than  this  dissertation’s  participation  factor.  To  find  the  relation¬ 
ship  between  the  two  find  the  displacement  of  the  sand,  mass  1,  and  fluid,  mass  2,  from  a 
first  mode  displacement  of  one: 
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Xii  =lxP, 


- 

X21  =  ^Pi 


(E-4) 


A 


The  displacement  of  the  fast  wave  from  a  unit  displacement  of  the  first  mode  is: 

/l=(l-/3>y„+/3Z,,  (E-5) 


Which  is  this  dissertation’s  participation  factor  eq.  (2-48). 

The  last  column  of  the  Modal  Computation  Table  in  Fig.  E.4  contains  what  this 
dissertation  calls  effective  forces  {Va  is  the  input  velocity). 

One  of  the  biggest  differences  in  the  mathematics  between  DDAM  and  Biot  The¬ 
ory  is  that  the  eigenvalues  for  DDAM  are  the  natural  frequencies  squared,  whereas  in  the 
Biot  Theory  the  eigenvalues  are  the  wave  numbers  squared. 

As  mentioned  earlier,  DDAM  had  a  big  influence  in  developing  the  lumped  pa¬ 
rameter  model. 
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SHOCK  DESIGN  OF  SHIPBOARD  EQUIPMENT 
PART  I  -  DYNAMIC  DESIGN- ANALYSIS  METHOD 


R.  O.  Belsheim  and  G.  J.  O'Hara 


INTRODUCTION 


During  recent  years  the  Bureau  of  Ships,  through  its  laboratories  and 
contractors,  has  carried  out  an  extensive  program  to  decrease  vulnerability 
'of  submarines  and  ships  to  underwater- explosion  attack.  This  vulnerability 
has  been  considered  in  terms  of  hull  and  equipment*  shock  damage,  each  of 
which  may  cause  ultimate  loss  of  the  ship  (1),  Perhaps  because  shock 
damage  is  not  as  definite  a  concept  as  hull  damage,  or  because  the  war 
damage  reports  did  not  show  it  to  be  as  damaging  (2),  the  experimental 
.research  effort  has  been  directed  more  toward  the  understanding  of  hull 
damage,  A  strong  exception  to  this  trend  has  been  the  design  of  nuclear 
power  plant  components,  which  today  are  stronger  on  the  average  than  other 
equipment  as  a  result  of  extensive  tests  and  resultant  design  procedures 
(3-9).  Another  notable  exception  is  the  Polaris  missile  installation,  which 
is  designed  to  protect  against  relatively  severe  shock  motions,  even  though 
the  missile  itself  is  very  fragile  (10).  Overall,  however,  recent  predictions 
of  damaging  radii  for  pressure  hull  damage  (11)  and  equipment  shock  damage 

(12)  suggest  that  equipment  is  relatively  more  prone  to  damage  in  modern 
submarines.  A  similar  conclusion  in  the  case  of  some  surface  ships,  i.e.  , 
destroyers,  is  found  from  study  of  results  of  recent  tests'  in  the  Pacific 

(13) .  Only  full-scale  tests  on  operational  ships  and  submarines  can  give 
a  definite  evaluation  of  vulnerability. 

Even  though  many  problems  remain,  a  considerable  amount  of  research 
has  been  done  and  the  authors  believe  information  is  available  so  that  better 
design  analyses  can  be  made.  The  purpose  of  this  report  is  to  review  this 
material  and  to  propose  a  dynamic  design-analysis  method  for  evaluating 
shock  resistance  of  shipboard  equipment  which  is  essentially  linear  and 
elastic  and  which  does  not  rest  on  noise  or  vibration  mounts  which  often  are 
non-linear  and/or  will  bottom.  All  of  the  required  information  is  not  yet 
available,  since  only  realistic  underwater  explosion  tests  can  give  correct 
inputs,  but  the  missing  information  can  be  obtained.  It  is  not  intended  that 
this  design  analysis- method  supersede  present  laboratory  shock-test  and 
evaluation  procedures,  which  will  always  be  desirable,  but  it  should  help 

^Equipment,  as  used  here,  includes  machinery,  electronic  equipment,  and 
other  items  which  are  attached  internally  to  the  ship  by  some  intervening 
structure;  it  does  not  include  items  subjected  directly  to  the  shock  wave  in 
the  water. 


Figure  E.3 

Author’s  Page  of  NAVSHIPS  250-423-30 
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Step  2  -  Calculate  the  set  of  influence  coefficients  for  these  points  and 
form  the  influence  coefficient  matrix. 

Step  3  -  Using  the  method  of  matrix  iteration,  find  the  mode  shapes  and 
natural  frequencies  of  the  first  few  modes,  say  two  or  three,  (There  is  a 
numerical  example  worked  out  in  Appendix  A  for  the  first  three  modes  of  a 
uniform  cantilever  beam.  )  Of  course  the  design  analyst  may  use  other  proce¬ 
dures  he  is  more  familiar  with.  For  example,  formulas  for  the  computation 
of  the  mode  shapes  and  natural  frequencies  of  a  two  degree  of  freedom  system 
are  given  in  Appendix  B. 

Step  4  -  For  each  separate  mode  "a"  complete  the  following  table: 


Modal  Computation  Table  (Mode  a) 


Mass 

Number 

Mass  =  Mj 

Mode  — 

Shape  ~ 

MiXia 

MiXil 

MiXiaPaVa^^a  * 

I 

Mj 

^la 

^l^la 

2 

M2 

5^2a 

^2^2a 

M2X2^P^V^(U^ 

- 

- 

- 

- 

- 

- 

n 

Mn 

^a 

MnXna 

MnXnaP  a  V  a^^a 

E 

MiXia 

MiXil 

Mi  Xi  aP  a^  a^a 

Step  5  -  Calculate  the  participation  factors  for  each  mode  "a" 


_  1  MiXia 

Pa  =  i (6) 

1 

Step  6  -  Calculate  the  "effective  mass"  in  each  mode  by  means  of  the 
formula 


r?  H^ia]  ^ 

Ma  =  Pa  2  MiXia  =  ^  - - - - .  (7) 

i  X 


Figure  E.4 

Modal  Analysis  of  NAVSHIPS  250-423-30 
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Appendix  F: 

PLANS  OF  TRAY  AND  LID 

The  following  AUTOCAD  plans  show  the  construction  of  the  foundation,  tray, 
and  lid  used  in  the  reflection  ratio  method.  Fig.  F.  1  shows  the  assembly  with  the  receiver 
and  transmitter  under  the  sand.  Fig.  F.2  is  the  tray  that  holds  the  sand.  Fig.  F.3  is  the  lid 
used  when  the  Ethafoam  is  on  the  top  surface  of  the  sand.  And  Figs.  F.4  and  F.5  show 
the  profile  and  plan  view  of  the  aluminum  foundation.  These  items  were  constructed  in 
the  machine  shop  at  the  University  of  Washington’s  Apphed  Physics  Laboratory. 


positioning  system  roller  ledge 
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Figure  F.  1 

Assembly  View  of  Equipment  in  Tank 
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Figure  F.2 
Plan  of  Tray 


Rod  <4) 
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Figure  F.3 
Plan  of  Lid 


inside  dimension 
MTUDLE  rnUNDATIDN 

Plan  View 


Figure  F.5 

Plan  View  of  Foundation 
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Appendix  G: 

PULSE  FORMATION 

A  pulse  was  needed  that  had  no  side  lobes  after  0.03  ms.  This  is  the  time  for 
sound  to  go  through  4”  of  aerylie.  The  desired  reeeived  pulse,  y2,  eq.  (G-1),  is  shown  in 
Fig.  G.l.  The  desired  pulse  is  eomprised  of  a  sine  wave  at  250,000  Hz  with  a  Gaussian 
envelope. 


The  equation  for  the  desired  reeeived  pulse  eentered  at  0.35  ms  is: 

y2=exp  (t-0.35-10~^)  (-3.6968el0)  sin(2.;r-2.5e5T) 


(G- 
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To  find  the  transfer  funetion  of  the  system  the  transdueers  were  exeited  with  a 
pulse  xi,  eentered  at  0.05  ms,  shorter  in  time  than 72,  henee  with  more  frequeney  eontent. 

xi  =exp  4-0.05-10“^f  (-3.6968el0) 


sin(2.;r-2.5e5-t) 


(G-2) 


The  reeeived  signal  from  an  exeitation  of  xi  is  yi.  The  signal  yi  is  a  signal  refieeted  from 
a  water-air  interfaee  with  a  travel  distanee  seleeted  so  that  the  pulse  returns  at  0.35  ms. 

The  goal  is  to  manufaeture  an  exeitation  signal  for  the  transfer  funetion  —  that  will  pro- 

XI 

duce  the  desired  reeeived  signal  72-  This  was  done  by  eonverting  xi,  yi,  and  72  into  the 
frequeney  domain  and  finding  the  frequeney  eontent  ofx2. 

Z2=iF2  (G-3) 

Where  X2  is  the  Fourier  transform  of  X2,  etc.  The  frequency  contents  of  xi  and  yi  are 
shown  in  Fig.  G.2. 

The  forcing  function  fed  into  the  transmitting  hydrophone,  X2,  to  produce  the  de¬ 
sired  received  pulse,  72,  is  the  inverse  Fourier  transform  of  X2.  See  Fig.  G.3. 

Using  the  excitation  signal  X2  the  received  pulse,  _y2rec,  is  close  to  the  desired  re¬ 
ceived  pulse,  yi.  See  Fig.  G.4.  The  received  pulse  is  the  coherent  average  of  100  sam¬ 
ples  reflected  off  a  water-air  interface. 


Figure  G.2 

Frequency  Contents  of  Received  and  Excitation  Signals 


Amplitude 
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Appendix  H: 
SAND  SIZE 


Sand  size  and  other  properties  are  listed  below.  A  thanks  is  given  to  Kevin 
Briggs  of  the  Naval  Researeh  Laboratory  (NRL)  for  providing  the  sand  and  the  following 
properties.  The  sand  was  obtained  from  Ward’s  Natural  Seience  Establishment,  Inc., 
Rochester,  NY,  and  has  a  product  name  of  “Sand;  Ottawa,  Sea.”  The  sand  properties 
measured  at  NRL  (from  sand  with  the  same  product  name)  are: 


Phi  size  at  percentage  levels: 

5  16  25  50  75  84  95 

1.60  1.82  1.92  2.14  2.43  2.62  3.02 


Percentages  of: 
Gravel  Sand  Silt 

0.00  99.68  0.25 


Clay. 

0.08 


Moment  measures: 

Mean  S.Dev.  Skew  KG 

1.84  0.53  2.91  93.83 


Post- analytical  weight:  129.66 

Note  that  phi  units  are  logarithmic  such  that  the  grain  size  in  mm,  d,  is  related  to 
the  size,  (p,  in  phi  units  through  the  expression: 


d  =  2 


(H-1) 
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